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Abstract: We give the first full next-to- leading order analytical results in Chiral Per- 
turbation Theory for the charged Kaon K — > 3n slope g and decay rates CP-violating 
asymmetries. We have included the dominant Final State Interactions at NLO analytically 
and discussed the importance of the unknown counterterms. We find that the uncertainty 
due to them is reasonable just for Age, i.e. the asymmetry in the K + — > ir + ir + Tr~ slope g, 
we get Age = —(2.4 ± 1.2) x 10 -5 . The rest of the asymmetries are very sensitive to the 
unknown counterterms, in particular, the decay rate asymmetries can change even sign. 
One can use this large sentivity to get valuable information on those counterterms and 
on Im Gs coupling -very important for the CP-violating parameter e' K — from the eventual 
measurement of these asymmetries. We also provide the one-loop 0(e 2 p 2 ) electroweak 
octet contributions for the neutral and charged Kaon K — > 3ir decays. 
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1. Introduction 

The decay of a Kaon into three pions has a long history. The first calculations were done us- 
ing current algebra methods or tree level Lagrangians, see [1] and references therein. Then 
using Chiral Perturbation Theory (CHPT) [2, 3] at tree level in [4]. Some introductory 
lectures on CHPT can be found in [5] and recent reviews in [6]. 

The one- loop calculation was done in [7, 8] and used in [9], unfortunately the analytical 
full results were not available. Recently, there has appeared the first full published result 
in [10]. 

CP-violating observables in K — > 3ir decays have also attracted a lot of work since 
long time ago [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22] and references therein. 

At next-to-leading order (NLO) there were no exact results available in CHPT so 
that the results presented in [16, 17, 18, 19] about the NLO corrections were based in 
assumptions about the behavior of those corrections and/or using model depending results 
in [16]. In [20, 21] there are partial results at NLO within the linear a-model. 

Recently, two experiments, namely, NA48 at CERN and KLOE at Frascati, have 
announced the possibility of measuring the asymmetry Age and Agjy with a sensitivity of 
the order of 10 -4 , i.e., two orders of magnitude better than at present [23], see for instance 
[24] and [25]. It is therefore mandatory to have these predictions at NLO in CHPT. The 
goal of this paper is to make such predictions. 

In particular, we have explicitly checked the one-loop results of [10], we also provide 
the complete one-loop calculation for the electroweak octet contribution up to 0(e 2 p 2 ) in 
CHPT for all the decays K — > 37r and finally, we estimate the dominant FSI for the charged 
Kaon K — ► Sir decays. We use all this to make the first full NLO in CHPT predictions 
for the charged Kaon K — > 3ir slope g and decay rates CP-violating asymmetries. We also 
present analytical results for all of our predictions. 

Notation and definitions of the asymmetries are in Section 2. In Section 3 we collect 
the inputs that we use for the weak counterterms in the leading and next-to-leading order 
weak chiral Lagrangians. In Section 4 we give the CHPT predictions at leading- and next- 
to-leading order for the decay rates and the slopes g, h and k. We discuss the results for 
the CP-violating asymmetries at leading order first in Section 5 and we discuss them at 
NLO in Section 6. Finally, we give the conclusions and make comparison with earlier work 
in Section 7. In Appendix A, the AS = 1 CHPT Lagrangian used at NLO can be found. 
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In Appendix B we give the notation we use for the K — > 3tt amplitudes and the new results 
at order e 2 p 2 . In Appendix C we give the analytic formulas needed for the slope g and the 
asymmetries Ag at LO and NLO and in Appendix D the relevant quantities to calculate 
the decay rates T and the CP-violating asymmetries in the decay rates Ar also at LO and 
NLO. In Appendix E we give the analytical results for the dominant -two-bubble- FSI 
contribution to the decays of charged Kaons and to the CP-violating asymmetries at NLO 
order, i.e. order p 6 . 

2. Notation and Definitions 

The lowest order SU(3) x SU(3) chiral Lagrangian describing |A5| = 1 transitions is 



(2) 

|AS|=1 



CFl e 2 G E tr (a 32 u^Qu) + CF 4 [G 8 tr (A 32 u^) + G" 8 tr (A 32 X+) 



+ G 27 fi> kl ti(A ijUf ,) tr(A fc ^) 



+ h.c. 



with 



C = -l9± Vud V,: s ~ -1.07 x 10- 6 GeV~ 2 . 



(2.1) 
(2.2) 



The correspondence with the couplings c 2 and c 3 of [7, 8] is 



C3 



c 2 — CFq Gs; 



(2.3) 



Fq is the chiral limit value of the pion decay constant f n = (92.4 ± 0.4) MeV, 

n M = iu ] {DJJ)v) = ul , 

Ajj = u\ijv) (Xij) a b = Siadjb , 

X+(_) = n f xn f + {-)ux ] u 



(2.4) 



X — diag(m ?t , m d ,m s ) a 3 x 3 matrix collecting the light quark masses, U = u 2 = 
exp {iy/2<&/Fo) is the exponential representation incorporating the octet of light pseudo- 
scalar mesons in the SU(3) matrix <5; 



/ id + IK 7j-+ 



$ = 



K+ 



\ 



TT 



V 



TT 1 V8 K 



The non-zero components of the SU(3) x SU(3) tensor t l ^ kl are 



^21,13 _ ^13,21 _ ^ . ^22, 23 _ ^23,22 



1 

"6 



^23,33 _ ^33, 23 _ _}_ . £23,11 _ ^11,23 _ 1 . 



(2.5) 
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and Q = diag(2/3, —1/3, —1/3) is a 3 x 3 matrix which collects the electric charge of the 
three light quark flavors. 

We calculate the amplitudes 

K 2 (k) - A^vrV) Ap 3 ) , [^oool. 

K 2 (k) - 7T + ( Pl )7T-(p 2 )7T°(p 3 ) , [A\_ Q ] , 

K^k) - 7r + ( Pl )7r-(p 2 )7r (p3) , [^+-o] > 
- ttV^VK+^s) , [^00+] , 
#+( fc ) 7r + ( Pl )7r + (p 2 )7r-(p 3 ) , , (2.6) 

as well as their CP-conjugated decays at NLO (i.e. order p 4 in this case) in the chiral 
expansion and in the isospin symmetry limit m u = rrid- We have also calculated the 
contribution of the 0(e 2 p 2 ) electroweak octet counterterms. In (2.6) we have indicated 
the four momentum carried by each particle and the symbol we will use for the amplitude. 
The states K\ and K 2 are defined as 

K m - ^ • (2.7) 

For the explicit form of the Lagrangian we have used, see Appendix A. Our results 
for the octet and 27-plet terms fully agree with the results found in [10] so that we do not 
write them again. The electroweak contributions to K — ► 3ir decays of order e 2 p° and e 2 p 2 
can be found in Subsection B.l in Appendix B. 

In this paper we discuss CP-violating asymmetries in the decay of the charged Kaon 
into three pions; namely, asymmetries in the slope g defined as 

jt^ 3 i Sl,g2,g3 l!2 = 1 + 9V + hy 2 + kx 2 + Oiyx 2 , y*) (2.8) 

and some asymmetries in the integrated K + — ► 3ir decay rates. Above, we used the Dalitz 
variables 

Sl- S 2 S 3 - S . . 

x = 7j — and y = s — (2.9) 

with Si = (k - p^ 2 , 3s = m 2 K + m^ (1) + m^ (2) + m^ (3) . 

The CP-violating asymmetries in the slope g are defined as 

g[K + — > 7r + 7r + 7r~] — g[K~ — > 7r~7r~7r + ] 



A 5 c = 



g[K + — > 7T+7T+7r J+^if — >• 7T 7T 7T+] 



, A S K+ ^ 7r°7r+ - g[K- 7r 7r vr- . . 

and A 57V = ^TTTT ^1T^ ^777 n — n T ■ 2 .10 

7T 7r vr+] + -> vr 7r°7r-] v ; 

A first update at LO of these asymmetries was already presented in [26]. 
The CP-violating asymmetries in the decay rates are defined as 

T[K + -► vr+vr+vr-] - F[K~ -► tt-tt--k + ] 



AT C = 



T[K+ -► vr+vr+vr-] + r[if- -► 71—71- Tr" 1 



r[K+ -► 7r 7r°7r+] - r[A"" -► 7r°7r°7r-] 
and A1 N = T[^+ - 7r07r07r+] + T[K- -+ ' (2 ' U) 
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In particular, we also want to check the statement that with appropriate cuts one can 
get one order of magnitude enhancement in AI^c and Ar^v asymmetries [13]. 



3. Numerical Inputs for the Weak Chiral Counterterms 

Here we collect the values of the weak chiral counterterms that we use in this work. 
3.1 Counterterms of the LO Weak Chiral Lagrangian 

In [10], a fit to all available K -> tttt amplitudes at NLO in CHPT [27] and K -> 3vr 
amplitudes and slopes in the K — > 3-ir amplitudes at NLO in CHPT was done. The result 
found there for the ratio of the isospin definite [0 and 2] K — ► tttt amplitudes to all orders 
in CHPT was 

A [K -> tttt] 



A 2 [K — »■ tttt' 

giving the infamous AI =1/2 rule for Kaons and 

A [K 7T7Tp (2) 



= 21.8; (3.1) 



A 2 [K -»■ 7T7T] 



17.8, (3.2) 



to lowest CHPT order p 2 . I.e., Final State Interactions and the rest of higher order cor- 
rections are responsible for 22% of the AI = 1/2 enhancement rule. Yet most of this 
enhancement appears at lowest CHPT order! The last result is equivalent [using Fq = 87.7 
MeV] to 

Re G 8 = 6.8 ± 0.6 and G 27 = 0.48 ± 0.06 . (3.3) 

In this normalization, Re Gs = G 2 j = 1 at large N c . No information can be obtained for 
Re (c 2 Ge) due to its tiny contribution to CP-conserving amplitudes. 

CP-conserving observables are fixed by physical meson masses, the pion decay coupling 
in the chiral limit Fq and the real part of the counterterms. To predict CP-violating 
asymmetries we also need the values of the imaginary part of these couplings. Let us see 
what we know about them. At large N c , all the contributions to Im Gs and Im (e 2 Ge) 
are factorizable and the scheme dependences are not under control. The unfactorizable 
topologies are not included at this order and they bring in unrelated dynamics with its 
new scale and scheme dependence, so that one cannot give an uncertainty to the large N c 
result for Im Gs and Im (e 2 Ge). We get 

ImG 8 =1.9 Imr, 

N c 

Im(e 2 G E ) =-2.9 Imr. (3.4) 

N c 

In the Standard Model [28] 

Im t = -Im ( Y?Y* j ~ -(6.05 ± 0.50) x 10~ 4 , (3.5) 
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and we used [29] 

(0|w|0) M g(2GeV) = -(0.018 ± 0.004) GeV 3 (3.6) 

which agrees with the most recent sum rule determinations of this condensate and of light 
quark masses -see [30] for instance- and the lattice light quark masses world average [31]. 

There have been recently advances on going beyond the leading order in l/N c in both 
couplings, Im Gg and Im (e 2 Gg). 

In [32, 33, 34], there are recent model independent calculations of Im (e 2 G^;). The 
results there are valid to all orders in 1/N C and NLO in as- They are obtained using the 
hadronic tau data collected by ALEPH [35] and OPAL [36] at LEP. The agreement is quite 
good between them and their results can be summarized in 

Im(e 2 G £ ) = -(4.0±0.9)Imr, (3.7) 

where the central value is an average and the error is the smallest one. In [37] it was used 
a Minimal Hadronic Approximation to large N c to calculate Im (e 2 Gs), they got 

Im (e 2 G E ) = -(6.7±2.0)Imr, (3.8) 

which is also in agreement though somewhat larger. There are also lattice results for 
Im (e 2 Gg) both using domain-wall fermions [38] and Wilson fermions [39]. All of them 
made the chiral limit extrapolations, their results are in agreement between themselves and 
their average gives 

Im (e 2 G E ) = -(3.2±0.3)Im r. (3.9) 

There are also results on Im Gs at NLO in 1/N C . In [40], the authors made a calculation 
using a hadronic model which reproduced the AI =1/2 rule for Kaons through a very large 
Q2 penguin-like contribution -see [41] for details. The results obtained there are 

Re G 8 = 6.0 ± 1.7, and G 27 = 0.35 ± 0.15 , (3.10) 

in very good agreement with the experimental results in (3.3). 
The result found in [40] is 

Im G 8 = (4.4±2.2)Im r (3.11) 

at NLO in l/N c . The hadronic model used there had however some drawbacks [42] which 
have been eliminated in the ladder resummation hadronic model in [43]. The work in 
[40, 41] will be eventually updated using this hadronic model. 

In [40] there was also a determination of Re (e 2 GE) though very uncertain. However, 
since the contribution of Re (e 2 Gg) is very small in all the quantities we calculate, we take 
the value from [40] with 100% uncertainty and add its contribution to the error of those 
quantities. 

Very recently, using a Minimal Hadronic Approximation to large N c , the authors of 
[44] found qualitatively similar results to those in [40]. I.e. enhancement toward the 
explanation of the AI = 1/2 rule through Q2 penguin- like diagrams and a matrix element 
of the gluonic penguin Qq around three times the factorisable contribution. The same type 
of enhancement though less moderate was already found in [45]. 
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3.2 Counterterms of the NLO Weak Chiral Lagrangian 

To describe K — ► 3ir at NLO, in addition to Re G 8 , G27, Re (e 2 G^), Im G 8 and Im (e 2 Gs), 
we also need several other ingredients. Namely, for the real part we need the chiral logs 
and the counterterms. The relevant counterterm combinations were called Ki in [10]. The 
chiral logs are fully analytically known [10] -we have confirmed them in the present work. 
The real part of the counterterms, Re Ki, can be obtained from the fit of the K — > 37r 
CP-conserving decays to data done in [10]. The relation of the Ki counterterms and those 
defined in Appendix A, and the values used for them are in Table 1 and Table 2 respectively. 





Re (G 8 )(7V£ - 27V? + 2N r 8 + JVJ) + G 27 {-\D$ 


K 2 


Re (G 8 )(N[ + iVf ) + G 27 (\D^ - |^ 8 ) 


K 3 


Re (G 8 )(iV£) + G 27 (|^ 7 + |^ 8 ) 


K A 


G 27 (D\ -Dl + W r 7 ) 


K 5 


G 27 (D r 30 + D r 31 +2D r 28 ) 


K e 


G 27 (8D^ 8 - D r 2% + D r 30 ) 


K 7 


G 27 (-4D r 28 + D r 29 ) 


K s 


Re (G 8 )(2iV£ + 4iVf + JVJ - 2iV[ - 4/^ - 2iV[ 2 ) + G 27 (-§D[ + 


Kg 


Re (G 8 )(iV- + AT- + JVJ) + G 27 (-l£>g) 


K w 


G 27 (21^ - 2D£ - DTj) 


K n 


G 27 D' 7 



Table 1: Relevant combinations of the octet N[ and 27-plet D\ weak counterterms for K — > 37r 
decays. 





Re i^(M p ) from [10] 


Im Ki(Mp) from (3.13) 


K 2 (M p ) 


0.35 ± 0.02 


[0.31 ±0.11] Im r 


K 3 (M p ) 


0.03 ±0.01 


[0.023 ±0.011] Im r 


K 5 (M p ) 


-(0.02 ±0.01) 





K 6 (M P ) 


-(0.08 ±0.05) 





K 7 (M p ) 


0.06 ± 0.02 






Table 2: Numerical inputs used for the weak counterterms of order p 4 . The values of Re Ki and 
Im Ki which do not appear are zero. For explanations, see the text. 

For the imaginary parts at NLO, we need Im Gg in addition to Im G 8 and Im (e 2 G^). 
To the best of our knowledge, there is just one calculation at NLO in 1/N C at present [40]. 
The results found there, using the same hadronic model discussed above, are 

ReG 8 = 0.9±0.1 and Im Gg = (1.0 ± 0.4) Im r . (3.12) 

The imaginary part of the order p 4 counterterms, Im Ki, is much more problematic. They 
cannot be obtained from data and there is no available NLO in 1/N C calculation for them. 
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One can use several approaches to get the order of magnitude and/or the signs of 
Im Ki. Among these approaches are factorization plus meson dominance [47]. If one 
uses factorization, one needs couplings of order p 6 from the strong chiral Lagrangian for 
some of the Ki counterterms, see also [48]. Not very much is known about these 0(p 6 ) 
couplings though. One can use Meson Dominance to saturate them but it is not clear that 
this procedure will be in general a good estimate. See for instance [49] for some detailed 
analysis of some order p 6 strong counterterms obtained at large N c using also short-distance 
QCD constraints and comparison with meson exchange saturation. See also [50] for a very 
recent estimate of some relevant order p 6 counterterms in the strong sector using Meson 
Dominance and factorization. 

Another more ambitious procedure to predict the necessary NLO weak counterterms is 
to combine short-distance QCD, large N c constraints plus other chiral constraints and some 
phenomenological inputs to construct the relevant AS = 1 Green functions, see [43, 49, 51]. 
This last program has not yet been used systematically to get all the AS = 1 counterterms 
at NLO. 

We will follow here more naive approaches that will be enough for our purpose of 
estimating the effect of the unknown counterterms. We can assume that the ratio of the 
real to the imaginary parts is dominated by the same strong dynamics at LO and NLO in 
CHPT, therefore 

Im Ki Im Gs, Im , , 

TT-Fr - ir~7^ - 0-9 ±0.3 Imr, 3.13 

Re Ki Re G 8 Re G' 8 

if we use (3.11) and (3.12). The results obtained under these assumptions for the imaginary 
part of the Ki counterterms are written in Table 2. In particular, we set to zero those Im Ki 
whose corresponding Re Ki are set also to zero in the fit to CP-conserving amplitudes done 
in [10]. Of course, the relation above can only be applied to those Ki couplings with non- 
vanishing imaginary part. Octet dominance to order p 4 is a further assumption implicit in 
(3.13). The second equality in (3.13) is well satisfied by the model calculation in (3.12). 

The values of Im Ki obtained using (3.13) will allow us to check the counterterm 
dependence of the CP-violating asymmetries. They will also provide us a good estimate of 
the counterterm contribution to the CP-violating asymmetries that we are studying. 

We can get a second piece of information from the variation of the amplitudes when 
Im Ki are put to zero and the remaining scale dependence is varied between M p and 1.5 
GeV. We use in this case the known scale dependence of Re Ki together with their absolute 
value at the scale v = M p from [10]. 



4. CP-Conserving Observables 

Here we give the results for the CP-conserving slopes gc, he, and kc and the decay rate 
Tc of K + — > 7r + 7r + 7r~ and slopes §n, ^at, and k^ and decay rate Tn of K + — ► 7r + 7r 7r 
within CHPT at LO and NLO. These results are not new -see [10] and references therein- 
but we want to give them again, first as a check of our analytical results and second, to 
recall the kind of corrections that one expects in the CP-conserving quantities from LO to 
NLO for the different observables. 
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We will use the values of Re G$ and G27 in (3.3), and disregard the EM corrections 
since we are in the isospin limit and they are much smaller than the octet and 27-plet 
contributions. For the real part of the NLO counterterms, we will use the results from a 
fit to data in [10]. So, really these are just checks. 

The values of the NLO counterterms given in [10] were fitted without including CP- 
violating contributions in the amplitudes, i.e., taking the coupling Gg and the counterterms 
themselves as real quantities. The inclusion of an imaginary part for these couplings does 
not affect significantly the CP conserving observables. 

To be consistent with the fitted values of the counterterms of the 0(p 4 ) Lagrangian 
we do not consider any 0(p 6 ) contribution to the amplitudes in this section. Indeed, these 
counterterms, fixed with the use of experimental data and order p 4 formulas, do contain 
the effects of higher order contributions. We also use the same conventions used in [10] for 
the pion masses, i.e., we use the average final state pion mass which for K + — ► 7r + 7r + 7r~ is 
m n = 139 MeV and for K + — > 7r°7r°7r + is m n = 137 MeV. In the following subsections we 
provide analytic formulas at LO and in Tables 3 and 4 we give the numerical results. 

4.1 Slope g 

The slope g is defined in equation (2.8). We give here the results for 

9C = ^ \g[K + -»■ 7T + 7T + 7r~] + g[K~ -> 7T~7r~7T + ] j 

and g N = \{g[K + -► ttW] + g[K~ -► ttVtT]}. (4.1) 

Without including the tiny CP-violating effects g[K + ]^o = g[K~\]_,o, 

lo o 2 3Re Gg — I3G27 

9c = - 3rri 



T m\ (3Re G 8 + 2G 27 ) + 9F 2 Re (e 2 G E ) ' 
L() ml (19m 2 K - Aml)G 27 + 6(m 2 K - m 2 n )Re G 8 + 9F 2 Re (e 2 G E ) 



9n 



(m 2 K - m 2 ) m\ (3Re G 8 + 2G 27 ) + 9F 2 Re {e 2 G E ) 



(4.2) 

The value for Re (e 2 G^) is not very well known. However its contribution turns out to 
be negligible and for numerical purposes we take the result for Re (e 2 G_s) from [40] with 
100% uncertainty. We do not consider its contribution for the central values in Table 3 
and we add its effect to the quoted error. In addition, the quoted uncertainty for g c ° and 
g^f contains the uncertainties from Re Gg and G27 in (3.3). 

The analytical NLO formulas are in (C.2). It is interesting to observe the impact of 
the counterterms so that we calculate also the slopes at NLO with Ki = 0, see Table 3. 
The contribution of the counterterms at /j, = M p is relatively small for go and g^, see 
Table 3. 

4.2 Slopes h and k 

We can also predict the slopes /^(N) an d &c(iv) defined in (2.8). At LO, the slope kc 
for K + — ► 7r + 7r + 7r~ and the slope A;^ for K + — ► 7r°7r°7r + are identically zero and the 
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yc 


no~ 18 GeV) 

-L fj I X\J V.J V V ( 


yN 


r,, fin- 18 GeV) 


T f"\ 


n i a _i_ n no 
— U.lo ± U.Uz 


l.z ± U.z 


U.00 ± U.U4 


U.o ( ± U.U 1 


irom laDie l 


n 99 -U n 09 
— U.ZZ zt U.UZ 


O.l It U.D 


n ri -i- n 
U.Ol It u.uo 


u.yo it u.zu 


NLO, 
Ki(M p ) = 


-0.28 ±0.03 


1.3 ±0.4 


0.80 ±0.05 


0.41 ±0.12 


PDG02 


-0.2154 ±0.0035 


2.97 ±0.02 


0.652 ±0.031 


0.92 ±0.02 


ISTRA+ 






0.627 ±0.011 




KLOE 






0.585 ±0.016 


0.95 ±0.01 



Table 3: CP conserving predictions for the slope g and the decay rates. The theoretical errors 
come from the variation in the inputs parameters discussed in Section 3. In the last three lines, we 
give the experimental 2002 world average from PDG [52], and the recent results from ISTRA+ [23] 
and the preliminary ones from KLOE [25] which arc not included in [52]. 



corresponding slopes hc(N) are equal to The NLO results are written in Table 4 

together with the slopes obtained when the counterterms Ki are switched off at fi = M p . 
We can see that the slopes hc(N) an d kc(N) are dominated by the counterterm contribution 
contrary to what happened with gc{N) which get the main contributions at LO. 





he 


kc 


h N 


k N 


LO 


0.006 ± 0.001 





0.075 ± 0.003 


1 


NLO, Ki(Mp) 
from Table 2 


0.012 ±0.005 


-0.0054 ± 0.0015 


0.069 ±0.018 


0.004 ± 0.002 


NLO, 

Ki(M„) = 


0.04 ± 0.01 


0.0004 ± 0.0025 


0.15 ±0.05 


0.008 ± 0.002 


PDG02 


0.012 ±0.008 


-0.0101 ±0.0034 


0.057 ±0.018 


0.0197 ±0.0054 


ISTRA± 






0.046 ± 0.013 


0.001 ± 0.002 


KLOE 






0.030 ±0.016 


0.0064 ± 0.0032 



Table 4: CP conserving predictions for the slopes h and k. The theoretical errors come from 
the variation in the inputs parameters discussed in Section 3. In the last three lines, we give the 
experimental 2002 world average from PDG [52] , and the recent results from ISTRA+ [23] and the 
preliminary ones from KLOE [25] which are not included in [52]. 



4.3 Decay Rates 

The decay rates K — > 3ir with two identical pions can be written as 

j r s Zmax r s lmax 

with 
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simvn = (E* + E{) 2 - LfEf - m) + \[e? - mfj , 

S3max = K - mf and s 3min = (m* + rrij) 2 . (4.4) 

The energies E* = (s 3 — m 2 + m 2 ) / and £J* = (m^ — S3 — mf)/(2y/~s^) are those of 

the pions 7T 7 and ir l in the S3 rest frame. It is useful to define 

\A C \ 2 = ^{\A (K + vr+vr+vr-) | 2 + \A (K~ 7^7^+) | 2 } , 

|^| 2 = - vrW) | 2 + |4 (K~ - 7r°7r°7r-) | 2 } . (4.5) 

At LO and again disregarding the tiny CP-violating effects we get 
1 4LO12 _ I 4 LO 12 _ I 4 lo 12 _ 

l A c I = \ A ++-\ — \ A — +1 — 



ICfx 



Re G 8 (s 3 - m\ - m 2 ) + ^ (l3m 2 + 3rr& - 13s 3 ) + Re (e 2 G B ) (-2F 2 ; 



l<T - l^o° + | 2 = l^o°-| 2 = Id 2 ReG 8 (m 2 -s 3 ) 

+ a( 2 27 2 n ( 5 "4 + 19m 2 m 2 , - 4m* + s 3 (4m 2 - 19m 2 ,)) 

+Re 2fn/-n^ ^ ~ m « ~ 3S3 ) " 

L \" L K '""ir) 



(4.6) 



The amplitudes |A?(iV)| 2 needed for the NLO prediction are in (D.6) in Appendix D. 

The results for Tc and Tj^ at LO and NLO are in Table 3. The contribution of 
Re (e 2 Gg) is very small (around 1%) and we include it in the final uncertainty as in 
Section 4.1 together with the rest of input uncertainties. We have also included in Table 
3 the results with the counterterms K \ = at fi = M p . We can conclude from them that 
the decay widths are strongly dependent on the NLO counterterms contribution. 

5. CP- Violating Predictions at Leading Order 

The numerators of the asymmetries in (2.10) and (2.11) are proportional to strong phases 
times the real part of the squared amplitudes. At LO in CHPT, the strong phases start 
at one-loop and are order p A jp 2 while the real part of the squared amplitudes are order 
(p 2 ) 2 . The denominators are proportional to the real part of the squared amplitudes which 
are order (p 2 ) 2 , so the asymmetries (2.10) and (2.11) for the slope g and decay rates T are 
order p 2 in CHPT. 

We have checked that the effect of Re (e 2 G^) is very small also for the Ag and Ar 
asymmetries. For the numerics, we have used Re (e 2 G^) = and used the value in [40] 
with 100% variation to estimate its contribution which we have added to the quoted final 
uncertainty of the asymmetries. For the Re Gs and G27 we have used always the values 
in (3.3). For Im Gs and Im (e 2 G^), we have used two sets of inputs; namely, the large N c 
limit predictions in (3.4) and the values in (3.11) and (3.7). For the pion masses we have 
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used the same convention used in [10] and given here in Section 4. The results are reported 
in Table 5. 



5.1 CP Violating Asymmetries in the Slope g 

At LO, the CP-violating asymmetries in the slope Ag c ^ can be written as [26] 

2 

A 9c?N) - S B C(N) Im G 8 + D C(N) Im (e 2 G E ) , (5.1) 

where the functions -Bc"(iV) an d ^c(at) on ly depend on Re Gs, G27, mx and m n and can 
be found in (C.4) and (C.5) in Appendix C. Numerically, 

Ag^° ~ [1.16 Im G 8 - 0.12 Im (e 2 G B )] x 10~ 2 , 
Ag^° ~ - [0.52 Im G 8 + 0.063 Im (e 2 G E )] x 10" 2 . (5.2) 

From (5.2) and the inputs discussed in Section 3.1 we conclude that the asymmetries 





A 5 LO(10-5) 


ArL°(io- 6 ) 


A^°(10- 5 ) 


ArLO(io- 6 ) 


(3.4) 
(3.11) and (3.7) 


-1.5 
-3.4 ± 2.1 


-0.2 
-0.6 ±0.4 


0.5 
1.2 ±0.8 


0.8 
2.0 ±1.3 



Table 5: CP-violating predictions at LO in the chiral expansion. The details of the calculation 
are in Section 5. The inputs used for Im G 8 and Im {c 2 Ge) are in the first column. The difference 
between Ag^P here and the one reported in [26] comes from updating the values of Re Gs and G27 
from [10]. The error in the first line is not reported for the reasons explained in Section 3. 



Agc(N) are poorly sensitive to Im (e 2 G E ). This fact makes an accurate enough measure- 
ment of these asymmetries very interesting to check if Im G% can be as large as predicted 
in [40, 44, 45]. It also makes these CP- violating asymmetries complementary to the direct 
CP-violating parameter e' K where there is a cancellation between the Im G% and Im (e 2 G E ) 
contributions. 



5.2 CP- Violating Asymmetries in the Decay Rates 

The observables we study here were defined in (2.11). We can write them again as follows 

rssmax d CSlmax fa A U |2 

c(N) ~ i:r : i:r ax i^w 

where the extremes of integration are in (4.4), the quantities |^4o(iV) | 2 were defined in (4.5) 
and A | Aq(n)\ 2 are defined by 

A\A C \ 2 = ^{|A (K + iT+iT+iT-) | 2 - \A (K~ -► vr-7r-vr + ) | 2 } , 

/S.\A N \ 2 = 1{|A (K + - 7r°7r°7r + ) | 2 - \A (K~ - nW') | 2 } . (5.4) 
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At LO we get, 



AULO |2 



ImG 8 {G 27 (^^-<^ 4) ) 

+ Re (^(^-BW 
+Im (e*G E ) {ReGs(B^-B^C^) 



R (2) r (4) 
n E °27 



o(2)^(4) 



)} 



(5.5) 



where we do not show explicitly the Sj dependence of the functions B\ and C^' nor the 



subscript C or N in and B^' for the sake of simplicity. The analytical expressions for 
the functions B^ and are reported in Appendix D. 

In (5.5), we have used consistently the LO result for the denominator of (5.3) though 
its value is very different from the experimental number, see Table 3. 
The numerics for the asymmetries in the decay rates are 



,(4) 



Ar c ° ~ [0.24 Im G 8 + 0.03 Im (e 2 G E )] x 10~ 3 , 
Ar^° ~ - [0.88 Im G 8 + 0.13 Im {e 2 G E )] x 10~ 3 . 



(5.6) 



The results using the two sets of inputs discussed in Section 3 for Im Gg and Im (e 2 Gg) are 
reported in Table 5. The asymmetries in the width are also poorly sensitive to Im (e 2 G E ) 
thus also their accurate measurement will provide important information on Im Gg. 

In [13], it was noticed that the asymmetry Arc increases if a cut on the energy of the 
pion with charge opposite that of the decaying Kaon is made. Afterward, the authors in 
[14] claimed that if this cut is made at S3 = 1.1 x 4m 2 , the asymmetry is enhanced by one 
order of magnitude. We checked that the decay rate asymmetry Ar^ at LO changes from 
its value in Table 5 to Arc* = —5.6 x 10~ 6 , i.e. one order of magnitude enhancement when 
we perform such a cut in the integration, in agreement with [14]. It remains to see if this 
enhancement persists at NLO and how feasible is to perform this cut at the experimental 
side. We will come back to this issue in the conclusions in Section 7. This enhancement 
does not occur for Ar^r. 



6. CP- Violating Predictions at Next-to-Leading Order 

At NLO one needs the real parts at order p 4 , i.e. at one-loop, for which we have the exact 
expression, see Appendix B. To make the full discussion about CP-violating asymmetries 
at NLO in CHPT we also need the FSI at order p 6 that would imply to calculate K — ► 3ir 
amplitudes at two-loops. However, one can use the optical theorem and the one-loop and 
tree-level tttt scattering and K — > 3tt results to get the imaginary part of the dominant 
two-bubble contributions. The results for these dominant two-bubble FSI are presented in 
the next subsection. 
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6.1 Final State Interactions at NLO 

Though the complete analytical FSI at NLO are unknown at present, one can do a very good 
job using the known results at order p 2 and order p 4 for tttt scattering and for K ^ Sir 
together with the optical theorem to get analytically the order p 6 imaginary parts that 
come from two-bubbles. These contributions are expected to be dominant to a very good 
accuracy. We are disregarding three-body re-scattering since they cannot be written as a 
bubble resummation. One can expect them to be rather small being suppressed by the 
available phase space [22]. 

Making use of the Dalitz variables defined in (2.9) the amplitudes in (2.6) [without 
isospin breaking terms] can be written as expansions in powers of x and y, 

A ++ _ = (-2a! + a 3 ) - (Pi - \pz + Vhs) y + 0{y 2 , x) , 

A 00+ = -(-2ai + a 3 ) - (~Pi + \pz + v^Ts) y + 0(y 2 ,x), 
A 2 + _ = ( ai + a 3 ) R - (ft + 3 ) R y + 0(y 2 , x) , 
A\_ = (ai + a 3 ) / - (Pi + Ps) 1 y + 0(y 2 , x) , 
A 2 000 = 3 (ai + a 3 ) R + 0(y 2 ,x), 

A 1 000 = S(ai + a 3 ) 1 + 0(y 2 ,x), (6.1) 

where the parameters ai, Pi and ji are functions of the pion and Kaon masses, Fq, the lowest 
order AS = 1 Lagrangian couplings Gs, G' 8 , Gn, Ge and the counterterms appearing at 
order p 4 , i.e., L^s, K^s. We do not add here EM corrections since we expect them to be 
small and of the same size of isospin breaking effects in quark masses which we have not 
considered. The O(e 2 p ) and 0(e 2 p 2 ) contributions can be found in Subsection B.l. 

In (6.1), superindices R and / mean that either the real part of the counterterms or 
their imaginary part appear, respectively. In the remainder, the superscript (H — 0) will 
refer to the amplitude A(K° — > 7r + 7r~7r ) = (A\_$ + A+_ )/y/2), that is proportional to 
the full couplings and not only to the real or the imaginary part of such couplings. 

If we do not consider FSI, the complex parameters af R , Pf R and 7j NR -with the 
superscript NR meaning that re-scattering effects have not been included- can be written at 
NLO in terms of the order p 2 andp 4 counterterms and the constants £>i,o(i) = b\ 0(1) +-^1 0(1) 
and defined in (B.3), (B.4) and (B.9). They can be obtained from Appendix D by 

expanding the corresponding functions Bi, Ci and Hi as in (B.9). We get 

«r = -gs\b^- ] + \ e * (4r 0) -<^) 

i=27,E 

4 e K„ ,,+ -° , -4r + - , )«<. 

1=1,11 

°s*- e c4 (<;+-> +2 B r o, K e « ++ - , +^a )<+ -° ) )^. 

i=27,E i=l,ll 

P- R = -Gs B^- ] - I E G* 0) + 4t + _) " 

i=27,E 
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ft™ 



1=1,11 



i=27,E 



+1 e (^r^-^r^- 2 ^ ^' 



i=i,n 



V37s NR 



1 E + - \ £ (^ (++ - } + ^ (0 ° +) ) .(6.2) 

i=27,E i=l,ll 



At LO the expressions above give 



«3 



LO 



iC 
iC 
iC 
iC 



m 



K 



+ G 



m 



27- 



A" 



27 



+ e^F 2 ' 



10m 2 , 

-Go 



^27" 



27 



e 2 G B -F 2 



m, 



2 1 



-G s mi - G 27 -^ 



-G 



5m: 



27 



G 



,2 

'K 

5ml 



18(m^- - m% 



(5m| - 14m 2 ) + e 2 G E F ( 



3m: 



2(m^ - m 2 , 



27 



4(m| : — m 2 ~ 



(3m| - 2m 2 ) + e 2 G E F 2 



3m 2 



A{m? K — m 2 ) 



(6.3) 



with the constant C defined in (2.2). 

The strong FSI mix the two final states with isospin 1 = 1 and leaves unmixed the 
isospin 1 = 2 state. The mixing in the isospin 1 = 1 decay amplitudes is taken into 
account by introducing the strong re-scattering 2x2 matrix R [17]. The amplitudes in 
(2.6) including the FSI effects can be written as follows at all orders, 



A 



(i) 



(i) 

oo+ / R 



(2) 



A 



(2) 

000 / R 

A {2) \u 



I + i] 



A 
A 

'a 

A 



(i) 



(i) 

00+ / N R 

(2) \ 
+-0 
(2) 

000 / N R 



(1 + iS 2 )A^l_\ NR , 



(6.4) 



with the matrices 



1/11 
3 I 1 -2 



1 

3 



1 
-3 1 



(6.5) 



projecting the final state with 1=1 into the symmetric-non-symmetric basis [17]. The 
subscript R (NR) means that the re-scattering effects have (not) been included. In these 
definitions the matrix R, S 2 and the amplitudes A^ depend on si, s 2 and S3. 
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NR 



Up to linear terms in y, equation (6.4) is equivalent to 

f-^ + \ a A =(I + ,M)(- ai + f a3 

73,R = (1 + ^2 ) 73,NR • 

(6.6) 

Here, the matrix K and 52 are functions of the meson masses and the pion decay coupling. 
At lowest order in the chiral counting they are given by 



mLO 1 \m\ - 9ml ( Wml + 2m 2 K ) \ 

rf2vr^ \ m K + 6m n y 3 m K D (m ^ +3m 2 ){m ^_ 9m 2 ) / 

and 

rLO _i 2 / ro%- - 9ro% 

52 = -96^f^V<T^| (6 - 8) 

in agreement with [22]. 

If we substitute the values of the masses and the coupling constant Fq, we get 

LO= / 0.136 \ <5 2 LO = -0.050. (6.9) 

^0.050 -0.143 J 

We have also obtained the phase 82 and two combinations of the M NLO matrix 
elements at NLO when including the dominant FSI from two-bubbles obtained as explained 
before. The determination of all the elements of M NLO would require the calculation of 
the FSI at NLO for all the amplitudes in (6.1) -we only have done the charged Kaon 
decays. The analytical expressions for these NLO quantities are given in Appendix E.5. 
Numerically, we get 

NLO 

= 1 + i 0.156, 



(-<*! + §a 3 ) R 



{-Pi + \fo 



R 



NR 

NLO 



(-ai + ia 3 )i 

l + i0.569, and <5^ LO = -0.104. (6.10) 



+ 3&)nr 

6.2 Results on the Asymmetries in the Slope g 

As we have seen in Section 5, the electroweak contribution to Ag at LO proportional to 
Im (e 2 G_e) is at most around 10% of the leading contribution proportional to Gs while 
Re (e 2 G^) generates a negligible contribution. We include in our results the NLO absorp- 
tive part of the electroweak amplitude which is proportional to Im (c 2 Ge)- The rest of the 
electroweak amplitude is just used in the estimate of the errors. 1 . 

x The expressions for the order e 2 p° and e 2 p 2 contributions to all the decay K — * 37r amplitudes are in 
Appendix B.l. 
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In order to study the NLO effects in §c(n) an d &-9c(N) > it is convenient to introduce 
\ A{K + ^ 3vr) |2 = A + + yA + + (x,y 2 ), 

\A(K-^3tt)\ 2 = Ac + yA y + 0(x,y 2 ), (6.11) 



so that 



g [K + (-) -► 3vr] 



a: 



+(-) 



,+(-) 

^0 



A 5 



(6.12) 



Notice that the numerator and denominator in (6.12) are not the same as the difference 
g[K + -► 3tt] - -► 3tt] and the sum g[K + — ► 37r] ± gfi'T — > 37r] respectively. At NLO, 

the sum A^Ag +t4q~t4~ does not contain the FSI at NLO since they are part of the order p 6 
contributions, i.e. of the next-to-next-to-leading order effects for the real parts. However, 
the difference Aq — Aq A~ is proportional to the imaginary part of the amplitudes, 
therefore to have it at NLO we must take into account the FSI phases, i.e. we need to 
include the FSI at NLO only in the imaginary part. 



The analytical expressions of the functions Aq^ ^ and Ay ( ' at NLO are collected for 
the charged and the neutral Kaon cases in Appendix C. From these expressions, we get 
the following numerical results 



Ag^ LO ~ 0.66 Im G 8 + 4.33Im K 2 - 18.11 Im K 3 - 0.07 Im (e 2 G E ) 



x 10" 



A 5 



NLO 
N 



0.04 Im G 8 + 3.69 Im K 2 + 26.29 Im K 3 + 0.05 Im (e 2 G E ) 



x 10 



-2 



(6.13) 



Where we have used the values for Re Ki from the fit to CP-conserving K — > 3tt amplitudes 
[10]. The NLO counterterms Im Gs, Im (e 2 G^) and Im K$ are scale independent. In 
(6.13), we have fixed the remaining scale dependence from Im K 2 at ji = M p . For the 
only two unknown counterterms Im K 2 and Im K3, we have made two estimates of their 
effects. First, using (3.13) as explained in Section 3. The other estimate of the effects of 
Im K 2 and Im K3 is to put them to zero and to vary their known scale dependence between 
|U = M p and (i = 1.5 GeV. We include the induced variation as a further uncertainty in 
our predictions. 

Our final results for the slope g asymmetries at NLO are in Table 6. 



The central 





A^LO (10 - 5) 


Ar£ LO (io- 6 ) 


A«# LO (10- 5 ) 


Ar^O( 10 -6) 


KiiMp) from Table 2 
Ki(M p ) = 


-2.4 ± 1.2 
-2.4 ± 1.3 


[-11,9] 
1.0 ±0.7 


1.1 ±0.7 
0.9 ±0.5 


[-9,11] 
4.0 ±3.2 



Table 6: CP-violating predictions for the slope g and the decay rates T at NLO in CHPT. The 
details of the calculation arc in Section 6. The inputs used for Im G$ and Im (e 2 GE) are in (3.11) 
and (3.7), respectively. 
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values are obtained with the input values in Table 2 and the uncertainty includes the 
uncertainties of Im G%, Im (e 2 G^), the uncertainties of the counterterms quoted in Table 
2, the variation due to the scale explained above and the error due to the electroweak 
corrections. 

The contribution of the order p 4 counterterms Im Ki to Age is around 25% using the 
values in Table 2 and the dominant contribution is the term proportional to Im G$. For 
Agw we find a much larger dependence on the values of the Im Ki. Of course, since Im Ki 
are unknown these results should be taken just as order of magnitude results, a factor of 
two or three could not be unreasonable for Age and Ag^. The contribution of Im (e 2 GE) 
is smaller than a 10% of the dominant one for both Age and Agjy. 

6.3 Results on the Asymmetries in the Decay Rates 

We also only include NLO absorptive electroweak effects proportional to Im (e 2 GE) for 
the same reasons explained in the previous subsection. The analytical functions |yl^^| 2 
and Aj^^^l 2 needed to obtain the asymmetries in (5.3) at NLO are given in (D.6). Also 
as explained in the previous subsection, one should consistently not include FSI at NLO, 
which are order p 6 , in the squared amplitudes j^^^l 2 since they are part of the next-to- 
next-to-leading order corrections. On the contrary, one has to include FSI at NLO in the 
differences A|yl^^| 2 since these differences are proportional to the FSI phases. 

The results obtained numerically from (D.6) in terms of the imaginary part of the 
counterterms are 



^NLO 



at n n lo ~ 



-2.8 Im G 8 + 49.2 Im K 2 + 103.6 Im K 3 + 0.2 Im (e 2 G E ) 
-3.11m G 8 + 45.7Im K 2 + 56.3lm K 3 + 0.12Im (e 2 G E ) 



x 10~ 3 . 



x 10~ 3 



(6.14) 



In both cases the final value of the asymmetry is strongly dependent on the exact value of 
the Im Ki due to large cancellations in the contribution proportional to Im Gg- If we use 
the uncertainties quoted in Table 6 for Im Ki, the induced errors in Arc and Ar^ are 
over 100%. In Table 6, we just quote therefore ranges for the two decay rates CP-violating 
asymmetries. 

7. Comparison with Earlier Work and Conclusions 
7.1 Comparison with Earlier Work 

The asymmetries Age and Ag^ have been discussed in the literature before finding con- 
flicting results. The rather large result 

\Ag c \ - |A#jv| ^ 140.0 x 10" 5 , (7.1) 



was found in [16]. 

The upper bounds 



\Ag c \ < 0.7 x 10- 5 , (7.2) 
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at LO and 

\Ag c \ < 4.5 x 10" 5 , (7.3) 

at NLO were found in [17]. The NLO bound was obtained making plausible assumptions 
since no full NLO result in CHPT was used. 
In [18], 

Ag c ~ -0.16 x 10" 5 , (7.4) 

was found at LO and 

Ag c ~ -(0.23 ± 0.06) x 10~ 5 and Ag N ~ (0.13 ± 0.04) x 10~ 5 (7.5) 

in [19] also at LO. The authors of [18, 19] also made some estimate of the NLO corrections 
and arrived to the conclusion that they could increase their LO result up to one order of 
magnitude. But again no full NLO calculation in CHPT was used. 

The asymmetries Arc and Ar n have also been discussed before and the results found 
were also in conflict among them: 

|Ar c | ~ 31.0 x 10~ 6 and IAIVI ~ 100.0 x 10~ 6 , (7.6) 

in [16] , 

Ar c ~ -0.04 x 10~ 6 , (7.7) 

in [18], 

Ar c ~ -(0.06 ± 0.02) x 10~ 6 and Ar^ ~ (0.24 ± 0.08) x 10~ 6 (7.8) 

in [19], and 

Ar c ~ -1.0 x 10~ 6 , (7.9) 

in [21] -where we have used sui(osm) — 0.85 [52]. The result in [16] was claimed to be 
at one-loop, however they did not use CHPT fully at one-loop. We find, in general, that 
the results in [16] are overestimated as already pointed out in [17, 18, 19, 20, 21]. See [17] 
where some explanations for this large discrepancy are discussed. 

The results in [17, 18, 19] were reviewed in [53]. They used factorizable values for 
Im Gs and Im Ge, i-e. the couplings in (3.4), therefore their results have to be compared 
with the first row in Table 5. The reason of the difference between their results and ours is 
due mainly to the fact that these authors obtain the value of Re Gs using the experimental 
value for the isospin 1=0 K — > irn amplitude Re ao- This amplitude Re ao contains large 
higher order in CHPT corrections. Corrections of similar size occur also in Im ao when 
considered at all orders. However the authors used analytic LO formulas for Im ao as well as 
for the imaginary parts of K — > 3ir amplitudes. This asymmetric procedure of considering 
the real parts of the amplitudes experimentally and the imaginary parts analytically just 
at LO leads to a value for Re Gs which is overestimated. Therefore the CP violating 
asymmetries at LO are underestimated in [17, 18, 19]. The same comments apply to the 
predictions of e' K in those references as emphasized in [48]. Our result in Tables 5 and 6 
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fulfill numerically the upper bound found in [17] for Age at NLO but not the upper bound 
found there at LO because of the same reason explained above. 

The results in [21] were obtained at NLO using the linear cr-model. Recently, there 
was an update of those results in [20]: 

Age ^ -(3.4 ±0.6) x 10 -5 , (7.10) 

at LO and 

Ag c ~ -(4.2 ±0.8) x 10~ 5 , (7.11) 

at NLO in the linear cr-model. It is , however, unclear from the text, the values used for the 
gluonic and the electroweak penguins matrix elements to get those results. Though the LO 
result in (7.10) agrees numerically with our result in Table 5, we do not agree analytically 
with the results in [20] when the author says that the electroweak penguins contribution at 
LO is as much as 34% of the gluonic penguins contribution. We find that the electroweak 
penguin contribution is one order of magnitude suppressed with respect to the gluonic one. 

7.2 Conclusions 

We have performed the first full analysis at NLO in CHPT of the CP-violating asymmetries 
in the slope g and the decay rate T for the disintegration of charged Kaons into three pions. 
We have done the full order p A calculation for K ^ 3tt and completely agree with the recent 
results in [10]. To give the CP-asymmetries at NLO, one needs the FSI phases at NLO 
also, i.e. at two loops. We have calculated the dominant two-bubble contributions using the 
optical theorem and the known one-loop and tree level results in Appendix E as explained 
in Section 6.1. Due to the small phase space available for the re-scattering effects of the 
final tree pions one expects the rest of the FSI to be very suppressed. We have included 
this contribution in our final numbers. As a byproduct, we have predicted the isospin 
1=2 FSI phase at NLO and two combinations of matrix elements of the isospin 1=1 FSI re- 
scattering matrix R at NLO. They can be found numerically in Section 6.1 and analytically 
in Appendix E.5. We have given analytical expressions for all the results in the Appendices 
B, C, D, and E. 

Our final results at LO can be found in Table 5 and at NLO in Table 6. If we use 
the counterterms in Table 2, we find NLO corrections of the expected size, i.e. around 
20%, for Age and Ag^- With those values for the NLO counterterms, the CP-violating 
asymmetry Age is dominated by the value of Im Gs while the rest of the CP-violating 
asymmetries studied here, namely, Ag^, ATc and Ar^r, are dominated by the value of 
Im K2 and Im K3. 

Of course, our results in Table 6 depend on the size of Im K2 and Im K3. If their 
values are within a factor two to three the ones in Table 2 then the central value of Age 
changes within the quoted uncertainties for it, while the central value for Ag^ doubles. 
The asymmetries in the decay rates Arc* and Ar n can change even sign if we vary Im K2 
and Im K3 within the uncertainties quoted in Table 2. Therefore, we have presented for 
them just ranges. 
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We partially disagree with references [17, 19, 53] when the authors claim that one 
could expect one order of magnitude enhancement at NLO in all the asymmetries studied 
here. We find that for Age and Ag^ the NLO corrections are of the order of 20% to 30% . 
Only Ar^ and Ar^r can vary of one order of magnitude and even change sign depending 
on the value of Im Ki and Im K3. We also find that Age can be as large as —4 x 10~ 5 
both at LO and NLO while in the conclusions of [17, 19, 53] it was claimed that any of 
these asymmetries could not exceed 10 -5 within the Standard Model. 

In Section 5.2, we found that making the cut proposed in [13, 14] for the energy 
of the pion with charge opposite to the decaying Kaon, there is one order of magnitude 
enhancement for Arc hi agreement with the claims in those references. This result is 
however valid for our LO calculation. It remains unclear whether the cut can provide a 
real advantage at NLO since in this case the cancellation among the various counterterm 
contributions can mask the effect. In addition, it remains to see how feasible is to perform 
this cut experimentally. We do not find this enhancement for Ar^r. 

The measurement of these CP- violating asymmetries by NA48 at CERN and/or by 
KLOE at Frascati and/or elsewhere at the level of 10~ 4 to 10 -5 will be extremely interesting 
for many reasons. The combined analysis of all four CP-violating asymmetries Age, Agjy, 
ATc and Ar^r can allow to obtain more information on the values of the presently poorly 
known Im G$, and the unknown Im K2 and Im K%. Due to the different dependence on 
these parameters, if the measurement is good enough, one can try to fix Im if 2 and Im K3 
from the measurement of the asymmetries Agjy, ATc and Ar^ which are dominated by 
the order p 4 counterterms and use them to predict more accurately Age- 

The large dependence of the asymmetry Age of Im Gs at NLO can also be used as 
consistency check between the theoretical predictions for Age and for the CP-violating 
parameter e' K . Any prediction for e' K has to be also able to predict the CP-violating 
asymmetries discussed here. In particular, the measurement of Age may also shed light on 
a possible large value for Im Gs as found in calculations at NLO in 1/N C -see for instance 
[40, 44, 45]. 

Moreover, it seems that some models beyond the Standard Model can reach values not 
much larger than 1 x 10 -4 for the CP-violating asymmetries, see for instance [54]. Our 
results can help to distinguish new physics effects from the Standard Model ones in these 
observables and unveil beyond the Standard Model physics. 
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Appendices 

A. AS = 1 Chiral Lagrangian 

At next to leading order, the SU(3) x SU(3) chiral Lagrangian describing K 
is given by 



3tt decays 



... A , ;1 = CF 2 Re G 8 {^Of + N 2 O s 2 + N 3 8 3 + iV 4 Of + N 5 Of + iV 6 Of + iV 7 Of 

+ N 8 Ol + N 9 Ol + N 10 Of + N u 8 n + N l2 0\ 2 + N 13 Of 3 } + h.c. (A.l) 

for the octet part [8, 55, 56], 



-(4) 



, : Vi| 1 = CF*G 27 {D^f + D 2 Of + D A Of + D 5 2 5 7 + D % Of + D 7 Of 

+ D 26 Of + D 27 Of 7 + D 2S Ofz + D^OllDsoOH + AnOf 7 } + h.c. (A.2) 

for the 27-plet part [8, 55] and 



-(4) 



|As| 1 = Ce 2 F 4 Re G 8 {Z 1 Of w + Z 2 0% w + Z 3 0§ w + Z 4 Of w + Z 5 Of w Z 6 Ol 

+ z 7 of w + z 8 oi w + z 9 o§ w + z w of w 

+ ZnO™ + Z l2 Of 2 w + Z 13 Of 3 w + Z ia O™} + h.c. 

for the electroweak part with the dominant octet structure [57]. 
The octet operators are 



EW 



(A.3) 



o 



O 



Of = tr ( AaatviA^iO 
Of = tr (A 32 n^^) tr (u^u") 
tr (A 32 (x+^ M + u^ X +)) 
Of = tr(A 32 x+)tr(nX) 
tr (A 32 (x-u^ - u^x-)) 
tr (A 32 x+)tr (x+) 
: tr(A 32 X-)tr (x_) 



n 8 



ro 8 

u 13 



Of = tr (A^u^uV) 



tr (A 32 u M ) tr (u u u^u v ) 



of 



0% = tr (A 32 ^) tr Kx+) 
Of = tr (A 32 «^) tr (x+) 
O 8 = tr(A 32 x+X+), 
= tr (A 3 2X-X-) , 



CO 8 
L/ 12 



(A.4) 



The 27-plet operators are 



27 



Of 

of 



of 
of 



f^tr(A lJ x+)tr(A MX+ ) 
f^tr(A, j x-)tr(A jfca _) 

t^tr (A^) tr (A kl (u» X+ + xW)) 
i^tr (A ijUfl ) tr (A fci («" x _ - x _u")) 
t^ fci tr(A iiX +)tr(A fc/ n^ At ) 
f^tr(A^)tr(A w ^)tr(x+) 
Of = i^ fc, tr (A^m) tr (Ahiu u u u ) 
Of = i«> fcl tr (A^ (^u„ + tr (A fc , + u u u»)) 



n 27 
u 6 



Of 
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O 27 



0% = t^tv (Ay«,) tr (A kl u v u»u v ) , 
t^ kl tr (A ijU ^ tr (A fc , {v?u v u v + , 

027 = t ^,fcl tr (A . tr (Afc ^) tr _ 



(A.5) 



The dominant octet electroweak operators are 

Of w = tr (A32 {u^Qu, x+}) , Of" = tr (A 32 u^Qu) tr ( X 

Of™' = tr (A 32 utQu) tr (x+^Qu) 

Of w = tr(A 32 «X) 
Of w = tr (A 32 «%) tr {QU^QU) 



O 



EW 



tr (A 32 n^) tr [Quu^ 



Of w = tr (A 32 {^Qu, ^}) tr (uQu^) 
Of 3 w = tr (A^tQn) tr (u%) 



0f w = tr(A 32X+ )tr {QrfQU), 
0£ w = ti (A 32 {^Qu,u^}), 
0§ w = tv(A 32 u»)ti (Qut^n), 

0f o w = tr (A 32 u") tr ({uQut, v)Qu] U(i ) , 

0f 2 w = tr (A 32 {u^Qu, u»}) tr (u^Quu^) , 

Of 4 w = tr (A 32 ntQn) tr (v)Quu^) . 

(A.6) 



We have done the NLO calculation in the presence of strong interactions which at LO 
order are described by [2, 3] 



= El [ tr {u ^) + tr ( X+ )] ; 



(A.7) 



at NLO, the SU(3) x SU(3) strong chiral Lagrangian needed in K — ► 37r decays is given 
by 

£( 4 ) = Litr {u^f + L 2 tr (it'V) tr (i^u") + L 3 tr {vPu^fuv) + L 4 tr (u^u^) tr ( X+ ) 

+ L 5 tr (u^x+) + L 6 tr ( X+ ) tr ( X+ ) + L 7 tr ( X _) tr ( X _) + ^L 8 tr (x+X+ + X-X-) • 

(A.8) 



B. K -> 3tt Amplitudes at NLO 

A general way of writing the decay amplitude for K + — ► 3tt at NLO including FSI effects 
also at NLO is 

A(K + -»■ 37r) (si,s 2 ,s 3 ) = G 8 a 8 (si,s 2 ,s 3 ) + G 27 a 27 (si, s 2 , s 3 ) + e 2 G £ a E (si, s 2 , s 3 ) 

+ F^(K i ,Z i ,s 1 ,s 2 ,s 3 ) + iF( 6 )(^,Z,, Sl , S2 , S3 ). (B.l) 

While for the corresponding CP conjugate the amplitude is 

A{K~ -> 37r) (si,s 2 ,s 3 ) = Gga 8 (si,s 2 ,s 3 ) + G 27 a 27 (si, s 2 , s 3 ) + e 2 G* E a E {s\, s 2 , s 3 ) 

+ F^\K*,Z*, Sl ,s 2 ,s 3 ) + iF^(K*,Z*, Sl ,s 2 ,s 3 ). (B.2) 
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The energies Sj are denned in Section 2, the Ki and Zi are counterterms appearing at 0(p 4 ) 
and 0(e 2 p 2 ) respectively, see Table 1 and (A. 3) for definitions. The functions F^ 4 ^(si, s 2 , s 3 ) 
and F( 6 )(si, S2, S3) are 

F^(k i ,Z t ,s 1 ,s 2 ,s 3 ) = H^( Sl ,s 2 ,s 3 )K t + J ( i 4) (s 1 ,s 2 ,s 3 )Z t , 

i=l,ll i=l,14 

F( 6 )(^,^, Sl , S2 ,s 3 ) = Hf\ Sl , S2 ,s 3 )K t + 4 6) {s 1 ,s 2 ,s 3 )Z i . (B.3) 

i=l,ll i=l,14 

The complex functions a% can be written in terms of real functions as 

ai(s 1 ,s 2 ,s 3 ) = Bi(si,s 2 ,s 3 ) + iCi(s!,S2,s 3 ) (B.4) 

Bi(s\, s 2 , s 3 ) and Cj(si, s 2 , s 3 ) are real functions corresponding to the dispersive and ab- 
sorptive amplitudes respectively and admit a CHPT expansion 

Bi( Sl ,s 2 , s 3 ) = B\ 2 \s u 8 2 , s 3 ) + Bf\s u s 2 , s 3 ) + 0(p 6 ) , 

Ci( Sl , S2 ,s 3 ) = C^( Sl ,S2, S3 ) + C^( Sl , S2 , S3 ) + 0(p 8 ), (B.5) 

where the superscript (2n) indicates that the function is 0(p 2n ) in CHPT. 

The functions B^ 2 \ b\ A \ cf^ and the part depending on Ki of in (B.l) and (B.5) 
for i = 8,27, which correspond to the CP-conserving amplitudes up to order 0(p 4 ) and 
without electroweak corrections, that is, 

A(K - 3vr) = Re G 8 (^ 2) ( Sl , s 2 , s 3 ) + B ( ^ ( Sl ,s 2 , s 3 ) + ^ 4) (s 1; s 2 , s 3 )) 

+G 2 7 (#27 ( S 1> S 2,S 3 ) + J B27 ) ( S 1' S 2,S 3 ) + iC^ ) (si,S 2 ,S 3 )) 

+FW(Re^,s 1 ,S2,s 3 ) (B.6) 

were obtained in [10]. We calculated these amplitudes for all the decays defined in (2.6) 
and got total agreement with [10]. The explicit expressions can be found there taking into 
account that the relation between the functions defined here and those used in [10] is, for 
the charged Kaon decays, 

M 10 (s 3 ) + M u { Sl ) + M n (s 2 ) + M 12 (si)(s 2 - s 3 ) + M 12 (s 2 )(s! - s 3 ) 
= ReG 8 (^ + ^+^ 4) ) (++ _ ) 

+ G 27 (< + < + ^?) (++ _ ) + ^ + „ ) , 
M 7 (s 3 ) + M 8 (si) + M 8 (s 2 ) + M 9 ( Sl )(s 2 - s 3 ) + M 9 (s 2 )(s! - s 3 ) 
= ReG 8 (^ + ^ ) + ^ 4) ) (oo+) 

+G 27 (*g> + < + *cf 7 } ) (Qo+) + F$ +) . (B.7) 

The functions B^\ B^\ and the part depending on Ki of F^ in (B.l) and (B.5) 
were calculated for i = 8, 27 in [10]. We calculated these quantities and got total agreement 
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with [10], the explicit expressions can be found there. The functions C^\si, s 2 , S3) (for 
i=8,27) and are associated to FSI at NLO coming from two loops diagrams and are 
discussed in Appendix E. 

We have also calculated the contributions of order e 2 p° and e 2 p 2 from the CHPT La- 
grangian in (2.1) and (A. 3) in presence of strong interactions for all the K ^ 3ir transitions, 
that fix the functions B^\ B^\ and jf^ . The results are in the next subsection. 

In order to calculate the asymmetries in the slope g defined in (2.8) we need to expand 
these amplitudes in powers of the Dalitz plots variables x and y, 

S1-S2 , S3 - S ,„ „s 

x = 5 — and y = * — . (B.o) 

The notation we are going to use here is 

G\ 2n \ Sl ,s 2 ,s 3 ) = G% n) + y G ff + 0(x,y 2 ) ; (B.9) 

where the function Gj(si, S2, S3) can be any of the functions £?i(si, S2, S3), Cj(si,S2,S3) 
defined in (B.4) or S2, S3), Ji(si,s 2 ,s 3 ) in (B.3). The coefficients G^ 2 ™^ are real 

quantities that depend on the masses m 2 , m 2 K , the pion decay constant and the strong 
counterterm couplings of 0(p 4 ), i.e., L\. 

B.l O(e 2 p ) and 0{e 2 p 2 ) Contributions 

Here we give the order e 2 p° and e 2 p 2 contributions to all the K —> 3ir amplitudes without 
including virtual photon loops ones. We do not give the order e 2 p 2 contributions from one 
order e 2 p° AS = vertex when it is inserted in the external Kaon and pion lines since these 
contributions are mainly responsible of the — m^ and m 2 K+ — m 2 K0 mass differences 
and we take physical masses for them. 

Electroweak interactions break in general isospin symmetry. An isospin decomposition 
which is valid up to order e 2 p 2 at least, is the following 

A + +-(si,s 2 ,s 3 ) = 2A Ci+ (si,s 2 ,s 3 ) + B c (s 1 ,s 2 ,s 3 ) + B t>+ (si, s 2 , s 3 ) , 

A 00+ (s 1 ,s 2 , s 3 ) = A C)0 (si, s 2 , s 3 ) - £ c (si, s 2 , s 3 ) + £ t ,+(si, s 2 , s 3 ) , 

2 

A\_ (s 1 ,s 2 , s 3 ) = Cq^Si, s 2 , s 3 ) + - [B^ (s 3 , s 2 , si) - B^ (s 3 , si, s 2 )] 

+A I n + (s 1 , s 2 , s 3 ) - B I n (s 1 ,s 2 , s 3 ) , 
2 

^+-o( s i) s 2, S3) = Cq(s 1 ,s 2 , S3) + - [£' (s 3 , S 2 , Si) - S/ (S3, Si, s 2 )] 

+A^ + (S1, S 2 , S3) - 5„ (Si, S 2 , S3) , 

Adoo( s i> s 2,S3) = 3A^ (si,s 2 ,s 3 ) , 

A 2 l00 (s 1 ,s 2 , s 3 ) = 3yl^ (si, s 2 , s 3 ) , (B.10) 

where A c , A n , B c and B n describe 1 = 1 final states, B t the 1 = 2 one and Co the 1 = 
one. In (B.10) the superindex R means that only the real part of the counterterms appear 
in these functions and the superindex / that only the imaginary part is present. In the 
following we give the order e 2 p 2 and e 2 p° contributions to the functions defined in (B.10). 
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In order to write the formulas in a concise form we define 

3Gp F 6 2 

A = m 2. - m £ , (B.fl) 
B£ p = B(m L ,m P , s k ) , £ LP = B(m L , m P , a) , (B.12) 
£>£ p = B 1 (m L ,m P ,s k ) , D tp = Ii(mi l m Pl s) , (B.13) 

where L, P = ir, K, r] and A; = 1, 2, 3. The functions B(m,L, mp, s k ) and B\(mi,, mp, s&) 
can be found, e.g., in [58]. The constants C ew must be used instead of C eu , in the functions 
with superindex i? or /, respectively. 

B.l.l Electroweak Contributions at O(e 2 p ) 

A n>+ (s 1 ,s 2 , s 3 ) = A n>0 (s 1 ,s 2 , s 3 ) = , 

A c ,+ (si,s 2 ,s 3 ) = A c fi(si,S2,s 3 ) = -C, 

C 

B n (s!,S 2 , S3) = + S 2 - 2s s) 1 

c 

B t ,+ (si,s 2 , s 3 ) = B t>0 (s 1 ,s 2 , s 3 ) = + S2- 2s 3 ) , 

5 c (si, s 2 , s 3 ) = + «2 - 2s 3 ) , 

C (si,s 2 ,s 3 ) = 0. (B.15) 

B.1.2 Electroweak Loop Contributions at 0(e 2 p 2 ) 

We checked that the electroweak loops at order e 2 p 2 do not break the isospin symmetry 
and 

A n ,+(s!, s 2 , s 3 ) = A nfi (si,s 2 , s 3 ) = A n (si,s 2 , s 3 ) , 
A c ,+(si, s 2 , s 3 ) = A cfi (s 1 ,s 2 , s 3 ) = A c (si,s 2 , s 3 ) , 
Bt,+ (si, s 2 , s 3 ) = B tj0 (si, s 2 , s 3 ) = B t (s u 32, s 3 ) • 



ew j 



(B.16) 



We get 



A n | Zj=0 = ^2 S (si) ,with 

■** 1=1,2,3 
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/A 2 

(sB KK + 4(m 2 - 8)3™) + — {m\ + 11m 2 ) - 7m% + 26m^ 



V s 



xm 2 + 29m^ - s(5m 2 K + 67m 2 ) + 27s 2 ) B Kw +8A (A + s) D Kw ] 



(B.17) 



A c | Zi=0 = ^{8[ml(-Ll + 8Ll + 4Ll) + ^Ll(m 2 K -3ml)] 



+ 



^ > with 

J=l,2,3 

Si(s) = [8(19^1 - 6m£)i/„ + 4(305777^ - 318m 2 ){v K + 

_ A(27A 2 + s(211m 2 K - 435m 2 ))(^ - v v ) - 6m 2 (-17m^ 

OS 

+5m 2 + % 8 )B™ + 162As5^ + 54(m^ - 10m 2 s + 3s 2 
+3m 2 Y (m 2 + s))B™ + 2(17m^ + 15m^ - m^(68m 2 - 43s) 



+ 17m 2 s - 24s 2 )B r > K + 



9A 3 



2(73m^ - 206m^m 2 + 25m^) 



+3(51m|- - 19m 2 )s + 16s 2 ) B Kn - 3(-9(5m|- - 6m|-m 2 
+ml) + A(m 2 K - 3m 2 )s + f>s 2 )D^ K + (261 x m A K + 153m^ 
-244m 2 s + 275s 2 - m 2 Y (414m 2 + 2\2s))D K ^ 



(B.18) 



Co\z i= o 
S 2 (si,s 2 ,s 3 ) 



f 2 

J 7T 

s_v 
216A 



S 2 (s 1 ,s 2 ,s 3 ) - S 2 (s 2 ,s 1 ,s 3 ) ,with 
6(2z^„ - i/jf - + 2f3A(si + s 2 - s 3 ) — — (3A 

V Si S9 S3 



+41s 3 )) K - - 9(s 3 - 4m 2 K )Bg K - 9(5mf Y - m 2 - s 3 ) 



x#J* + _( A +s 3 )(A-2s 3 ) J Bf " - 6AD r ^ K + 6(A - 253)!) 

S3 

+ 2lk i 9(4m ^ " Sl)(S2 " S3)B ^ " ( 42 "^ " 16m ' + ^ 
x(118m 2 - 81s 3 - 61si) + m 2 (9s 3 - 53«i) + 15si(3s 3 + si)) 

xBf + + l) (3A(si + s 2 - s 3 ) - (8ml- + 22m 2 - 3s 3 )si 

+Us 2 )B? n -6(m% - 3mi + 2m£(s 3 - 5si) + 2m^(m 2 - s 3 
-2si) + si(9s 3 + 2s 1 ))D n 1 K + 2(3A(si + s 2 - s 3 ) - (8m 2 K 



+22m 2 - 3s 3 )s! + Usl)D( 



Ktt 



(B.19) 



t\Zi=0 



ft lA 



l r 



-^3(53(51 + s 2 ) - 2s 1 s 2 ) + 2L r 5 (s 2 + si - 2s 3 )m 2 
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+ (S 3 (s!,s 2 ,s 3 ) + S 4 (s 1 ,s 2 ,s 3 ) +54(32,51,53)) I ,with 
Ss(s!, 8 2 , s 3 ) = —— IjT^&mjt + 39ml + 4m| f (7m 2 - 3s 3 ) - 30m 2 s 3 



1 


r ~ 9 


432A 


L(16tt 2 ) 


-3*1- 


6sis 2 ) 



2 / 108 

-— ( A 2 m 2 - 31m% - 26m^m 2 + 489m* + 3(49m|- 

A V S3 









^)) 







-(54(si + s 2 - 2s 3 )(2m 2 - sg)^ + 2(17rr& + 15m£ 
-68m 2 m|- + (17m 2 + 43m^)s 3 - 24sl)S^ - (— A 2 m 2 

+(61m^ + U2m 2 K ml + 13m 4 w ) - 6(31m 2 K + 51m 2 )s 3 + 185s 2 ;) 

x£f 71 + 3(45m 4 K + 9m£ - 54m|-m 2 + 4(3m 2 - m 2 K )s 3 -5a§) 
xZ)f + (189ml: + 81m^ - 270m^m^ 

+(404m 2 - 284m^)s 3 + 35s§)Df , (B.20) 



5 4 (5i, s 2 , 53) = ^ [54(9?r4 + m 2 , (3m 2 - Sl ) + ( Sl - 4m 2 )(s 3 + 2 Sl ))B^ 

+(Um% + 15ml + m 2 K (\3m\ - 54s 3 - llsi) - 64m 2 si + 3s x 
/ A 2 

x(18s 3 + Si))B? K - (9 — {m\ + 9m 2 - 2s 3 ) - m^ + 17m^m 2 
V si 

-16m* + 9{5m 2 K + m 2 )s 3 - 3(25m|- + 54m 2 - 3s 3 )si + 97s 2 ) 

x£f * - (-9A(8m 2 K - s 3 ) + 3(17m^ + 33m 2 - 27s 3 )si - 33s?) 
xDf - (-9A(10m|- - 6m 2 + s 3 ) + 5(23m|- - 53m 2 + 9s 3 )si 
+5s 2 )Z)H , (B.21) 



B, 



c\Z x =0 



( H [ ~ L ^ s 3(«i + to) ~ 2s lS2 ) + ^(8AL^ - 3m 2 Z7 5 ) 



f 2 

J TV 

x(si + s 2 - 2s 3 ) + (S 5 (s 1 ,s 2 ,s 3 ) + 5 6 (S1,S3) + S 6 (s 2 ,S3))} ' with 

55(51,52, s 3 ) = ^r^- 3^2 ( 5m l + 39m £ + 28m|-m 2 - 6s 3 (2m|- + 5m 2 ) 

-3{s 2 3 + 2sis 2 )) + 2(55m^ + 189(m 2 - s 3 ))u v + 2(175m| 

+513(m 2 - s 3 )){u K + u n ) + (3A(m^ - 7m 2 - 6s 3 )f- + — ) 

V \Si s 2 / 

+— (m|- + 8m 2 ) + -L(163mk + 554m^m 2 + 579m* - 243 
s 3 3A 

x (m\ + 3m 2 )s 3 )) (y K - v v ) + 4m 2 {-Ylm\ + 5m 2 + 9s 3 )Bf 
-108s 3 A#f K - 9(34m* - 42m^m 2 + (39m 2 K - 37m 2 )s 3 + 3s|) 
x5 3 ™ - \{l7m% + 15m* - GSm^m 2 + (43m^ + 17m 2 )s 3 - 24 
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xsl)B'l K + -(673m 4 K - 794m 2 Y m 2 + 337m^ - — A 2 (m 2 K 
6 S3 

+8ml) - 6(81m 2 K + 25m 2 )s 3 + 101s§)Sf 77 - ^(9(5m^ 
-Gm^m 2 + m*) - 4(rr& - 3m 2 )s 3 - 5s 2 )Df - ^(477m% 



+369ml + 1556m 2 s 3 - 445s 2 ; - m 2 K {UQml + 284s 3 ))-D : 



Ktt 



(B.22) 



5 6 (si,s 3 ) = 2^[2m^(17m|--5m2-9si)B? ,, -18(4ml- + ai(2s 3 + si)+4 
xm^m 2 - 2(s 3 + si)))#f K - 9(m^ + 3m^(9m 2 - 7«i) + 3s 3 si 

-m 2 (12s 3 - 851))^ + ^(332771^ - 177m£ + m 2 K {781ml - 594s 3 

1 A 2 

-311si) + 4^2(2753 - 13si) + 3(54s 3 - ll Sl ) Sl )Bf - -{9—{m 2 K 
-7ml - 6s 3 ) + (335tt4 + 104m^ - 27m 2 s 3 - m 2 K (439m 2 - 81s 3 )) 
-3(69777^ + 34m 2 - 27s 3 )si + 61s 2 )#f w + ^(9A(4(m|- + m 2 ) 

-3s 3 ) + (7777^ + 39m 2 + 27s 3 ) Sl - 49s 2 )D^ + i(9A(22777 2 Y - 34m 2 
+9s 3 ) - (U9m 2 K - 865m 2 + 189s 3 )si - 257s 2 )L>f H , (B.23) 



d 1 C ew c 2 

-Dn|Z i= 



£ 3 ( _s 3(si + s 2 ) + 2sis 2 ) - 2Lg(s 2 + si - 2s 3 )m 2 



ft l A 

+ (5 7 (si,S 2 ,S 3 ) + 58(51,53) + 58(52, S3)) } , with 

SV(si, 8 2 , s 3 ) = —fr-r T7r^(4m A K + Urn^-m 2 + 6m^ - 6s 3 (2m 2 K + m 2 n ] 



1 


- 6 


144A 


Ll6vr 2 


3(2s 2 


-si- 



-4(A(m 2 K - 7ml) (— - — - -\ - ^-(m 4 K + 5m 2 K m 2 n 

\ \S 3 Sl S 2 7 ^ 

+6m^ - s 3 (4m| + 5m 2 ))) - + - (l8(2s 3 - 8l - s 2 ) 



x (s 3 5 3 - (2m 2 + 83)5^) + (37m^ + 25m* - 134mj Y m 2 

? 3 J 



+ (54m 2 + 90m^)s 3 - 635^)5^+ 6((-19m 4 K + 2m 2 Y m 2 



+29m^) + — (m% - 9m^m 2 + 15m^m^ - 7m 6 w ) + (25m 2 K 
S3 

-19m 2 )5 3 - 3a§)Bf w + (135m^ + 6m| Y (-27m 2 + 53) + 3(9 
xm 2 - 5s 3 )(m 2 r + 3s 3 ))r»^ + 3(79m^ + 43m^ - 122m 2 m\ 
+ (110m 2 - 86m^)s 3 + 15s§)I)f , (B.24) 
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Ss(si, as) = i - 18 ( 2m K + "4( 6m - " 4s 3 - 3si) - Sl (3m 2 - s 3 - 2ai)) 

xSf x -18(m^(7m2 - Sl ) + 21m£ - 12m 2 (a 3 + si) + 3s 3 si)-Br 
+(-50tt4 + m^(-23m2 + 63s 3 ) + (m£ + 3si)(m 2 -9a 3 + 6«i)) 
xB^- 3(- 19m^ + 20m^- m|-(m2 - 6s 3 ) + — (m^ - 9m|m^ 

\ Si 

+\hm\m% - 7ml) + (^m 2 K + 5m 2 )si - 18a?) Bf w - 3(A(38 
xm^ - 26m 2 + 3s 3 ) + (-37m 2 K + 67m 2 + 3a 3 )ai - 21a?)L>f * 
-3(3A(8rr& - a 3 ) + (-5rr& - m 2 + 9s 3 ) Sl - 3s 2 )Z^] , (B.25) 

B.1.3 Electroweak Counterterm Contributions at 0(e 2 p 2 ) 

The electroweak counterterm contributions at 0(e 2 p 2 ) break isospin symmetry and we use 
the decomposition in (B.10). We define the constant 

Cg, = -4%KdCi^L Re G 8 . (B.26) 



'5^2 ua us f K fi 



We get 



A z fi ( Sl ,s 2 ,s 3 ) = l-2(m 2 K + 5m 2 7r )Z r 1 -2(2m 2 K + 5m 2 7r )Z r 2 + ^(m 2 K -2m 2 7r )Z r 3 
+ ™k ( _ 32 r _ 2 ^ + 6(Z[ 1 + Z[ 2 )) - (m 2 K - 27ml) f 
~(m 2 K - 3ml)(6Z[ 3 + Z r u )\ , (B.27) 



C z 2m 2 

A Ct+ ( Sl ,s 2 , a 3 ) = A C)0 (a!, a 2 , a 3 ) - -f ^ (-621 + 2Zf + Z r n + Z[ 2 ) , (B.28) 

5f( Sl ,a 2 ,a3) = ^(£l±^f^{-18(m^ + Ti4)^-18(2m^ + mi)^ 

+A(t2(Z 3 r - Z£) - 31Z 6 r + AZ r 7 + 9Z 8 r - 48Z[ 3 - 8Z[ 4 ) 

-9(3m|- - 2m\)Zl - Qm 2 K Z{ Q + 6(4m^ - 3m\)Z\ 1 

-6(4m 2 K - 5ml)Z[ 2 } , (B.29) 

S§,(»i,*2,*s) = + 1 S 2 2 A 253 {6(m 2 K + m 2 )Z[ + 6(2m 2 K + m 2 )^ 

-A(4Z£ + 3Z 6 r ) + 3(3m^ - 2m 2 )Z£ - 3(m 2 K - 2m 2 )Z£ 
+2m 2 K Z r w - 2mlZ[ 1 + 2(4m 2 K - 5m 2 )Z[ 2 } , (B.30) 
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Jit 



C% w Sl + s 2 - 2s 3 



(2^-4^ + ^), (B.31) 



A%J Sl ,S2,S 3 ) = 



r 1 



m 



K 



(3Z£ - 2Z r 6 + 2Z r 7 - 3Z r s - 3ZI - 2Z 



ft 9 

+2Z r u + 2Z{ 2 ) , 



(B.32) 



Af + (si,s 2 ,s 3 ) = A% ( Sl ,S2,S 3 ) + 



Cew Zmjc 
ft 3 



{-6ZZ + 2Z r 7 + Z\ x + Z[ 2 ) , (B.33) 



BZ( Sl ,s 2 ,s 3 ) = -^ff Sl+ H^ {lKrnl + ml)Z[ + l%{2rnl + ml)Zl 

-A(12Z 3 r + 6Z r 5 + 5Z 6 r - 8Z r 7 ) - 3(m 2 K - 4ml)(Z r 8 + Z r 9 ) 
-2(m\ - 4ml)Z r w + 2{\m\ - 7m^)(Z[ 1 + Z{ 2 )} . 



(B.34) 



C. The Slope g and A# at LO and NLO 

We have checked that the following relations 

• Fq Re (e 2 G E ) « m^Re G 8 

• Im G 8 « Re G 8 

• Im (e 2 G E ) « Re G 8 

can be used in this and the next sections to simplify the analytical expressions. To obtain 
the numerical results included in the text we use the full expressions, with no simplifications. 
We have also checked that the terms disregarded with the application of these relations 
generate very small changes in the numbers. 

Using the simplifications above, the value of g at LO can be written trivially as 



B^Re G 8 + B^Gm + B (2 [Re (e 2 G E ) 



>(2) 



(2) 



fl$Re G * + B 27,oG27 + B% Re (e^G E ) 



,(2) 



(2) 



(C.l) 



(2) (2) (2) 

The expressions for B% '■ , B 27 \ , and B y E ^ needed above can be obtained from the expressions 
of the corresponding B's for the charged Kaon decays + H — and 00+ in Appendix D and 
expanding them as in (B.9). The results we get are in (4.2). 

We consider now the NLO corrections to the slope g. Disregarding the tiny CP- 
violating we have g[K + — > 3ir] = g[K~ — ► 3ir] at NLO we get 



Aq LO 



E 

1=8,27 



8) 



Re Gi + 



3=1,11 



l j,o ' 



+ 



i=8,27 



C^Re G t ) 
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4 L ° = 2 



x 



+ 



E Re G i + £ fl$Re£ 

i=8,27 j=l,ll 
=8,27 j=l,ll 

2 (c*> <*) x £ (cg>Re Gi ) } 

=8,27 |_i=8,27 J J 



(C.2) 



for the coefficients defined in (6.11). 

One can get gc(N) a t NLO using (6.12) and the results above substituting the coeffi- 
cients cfj , Bf^ and H^f- by their values calculated expanding in the Dalitz variables the 
results in [10]. 

The slope g asymmetry in (2.10) can be written at LO as 



A 9c(N) ^ -£ B C (N) Im G 8 + D C(N) Im (e 2 G E ) ■ 



(C.3) 



We get 



Bq — -777— G271 

64 7T 



/ mt- — 9m 1 



m 



K 



+ 3ml 

.2 _2 



DC = 777 



\\m\m\ — 18m* + 5m^ 
X - m2)(3Re G 8 + 2G 27 )(13G 27 - 3Re G 8 ) ' 

3 1 / m 2 K — 9ml 



Mv\jm 2 K + 3ml m^(m| - m2)(3Re G 8 + 2G 27 )(13G 2 7 - 3Re G 8 ) 

[3Re G 8 (16m 2 K ml - 18m^ + 3m 4 K ) - G 27 (178m 2 K ml - 234ml + 69m^)] , 

(C.4) 



x 



and, in the neutral case, 



I m\ — 9ml 7m 2 + 4ml 



K 



m\ + 3ml 



E 



D N = 



9 1 \m 2 K - 9ml ml(18ml - 7m 2 K ] 



m 



K 



+ 3ml 



m K E 



+ H G lm K 



32 7T 



^271 



■? K — 9ml 3Qm* — \\9m\m\ — 60m^- 



m 



K 



+ 3ml 



with 



)Re G 8 ) 



(C5) 



(C6) 



E = (3Re G 8 + 2G 27 ) ((19m 2 K - 4m\)G 2 i + 6(m 2 K 

The sum Ay Aq + ^4^"^", necessary to get Ag at NLO using (6.11) and (6.12), can be 
obtained directly from (C.2) where we have neglected the small CP-violating effects. 
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For the difference A+A Q — Aq A y , we get 

{A^A, - A+A y ) NLO = AAi [(A 2 R - C 2 R ) V R - 2ArB r C r ] + 4Cj [{A 2 R - C 2 R ) B R 

+2ArPrDr] + 4 (B T C R - ViAr) {A% + C 2 R ) , (C.7) 

where A R , B R , C R and V R contain the contributions from the real parts of the counterterms 
= E (4? + 4?) ReG t+ Yl H S Re ^ . 

i=8,27,E i=l,ll 

** = E (4? + 4 4 i } ) Re ^ + E *£? Re * > 

i=8,27,E i=l,ll 
i=8,27,E ' i=l,ll 

^ = E (<$ + Re + £ R tl Re ^ • (C8) 

i=8,27,E i=l,ll 

While ^4/, £>/, C/ are the same expressions but substituting the real parts of the countert- 
erms by their imaginary parts. 

The coefficients B^ 2 ™^, and H^ 2 ™^ defined in (B.9) are real. 

D. The Quantities \A\ 2 and A\A\ 2 at LO and NLO 

Here we give the results for the quantities A and AA defined in (4.5) and (5.4), respectively. 
They enter in the integrands of the decay rates T in (4.4) and the CP-violating asymmetries 
Ar, see (5.3). 

To simplify the analytical expressions, we have made use of the fact that the imaginary 
part of the counterterms is much smaller than their real parts. The | A^^^ | 2 which give 
the asymmetries Ar at LO are in (4.6). 

The result for A|^4c| 2 at LO can be obtained substituting in (5.5) the functions 
Bf\s 1 ,s 2 , s 3 ) and G^(si, s 2 , s 3 ) for i = 8 and i = 27 by 

E>( 2 ) 1 \ ■ C Fo 1 2 2\ 

B y 8 l + _(s l ,s 2 , s 3 ) = 1 -3— (s 3 -m K - mj , 

Jir JK 

CK 4 1 



B 27++-( s i' s 2, s 3 ) = i : rff- ^(!3m 2 + 3m 2 K - 13s 3 ) , 

r W , > .CF 4 / 1 \ 

C8++ - (Sl ' S2 ' S3) = Z ^l-16^J 



x • 



- [s\ - s 3 {3ml + m 2 K ) + 2m 2 n (m 2 K + ml)} a(s 3 ) 



+ ^ [4s^ + s 2 (-4m 2 n + 2m 2 K - s 3 ) + m 2 n (4s 3 - 2m 2 K - 3m 2 .)] a(s 2 ) 
+ (exchange si and S2 in the second term) > , 
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n W , , -CF 4 f 1 



X • 



- [-13s§ + 3s 3 (13m2 + m 2.) - 2m^ (3m|- + 13m^)] <r(s 3 ) 



sl(14m* + 31tt&) + s 2 (26s 3 (m 2 K - ml) 



1 

36(772^ — m 2 ) 
-174m^m^ - 7m^ + 76m 4 ) + (104m£s 3 (m£ - m 2 K ) 
-168ml + 161m 2 K ml + 67m 4 K ml) a(s 2 ) 

+ (exchange si and S2 in the second term) > , (D.l) 



and 4 2) , 4 4) by 



R (2) _ ,- CF ( 9p 2\ 

BE++ --'m { 2F o ] 



+7T^ 2T t 3s 2 + s 2(5m^ - 12m 2 ) + m 2 (5m 2 - m 2 K )] a(s 2 ) 

4(m K — m%) 

+ (exchange s\ and S2 in the second term)| . (D.2) 

One can get A|Aat| 2 at LO substituting in (5.5) the functions Bf'\s\, s 2 , s 3 ) and 
Cf\si, s 2 , as) for i = 8 and i = 27 by 



B 8 00+( s l> S 2, S3) = i 73-p- (m 2 - S3) , 



C F 4 1 

5 27oo+( s i' S2 > s s) = i 73-7^ tt-^ 2T [5m 4 K + \§m\m 2 K - 4m 4 . + s 3 (4m 2 - 19m^)] , 

JttJk v(m K m^j 

^(4) , , . CF 4 / 1 \ 



X • 



- [s 2 3 + s 3 (m 2 K - ml) - mlm 2 K ] a(s 3 ) 



+ - [2s 2 2 + s 2 (s 3 - 2(4m 2 + m 2 K )) + m 2 (-4s 3 + 5m| + 9m 2 )] a{s 2 ) 
+ (exchange si and S2 in the second term) > , 



C 2700+(S1 , S2)S3 ) = %w - ^— ^ 



x • 



-1 
12 L 



26s 2 (m|- — m 2 ) + s 3 (56m 4 — 57m|-m 2 — \4m 4 K ) 
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+m 7T (31m K m n — 30m n + l9m K ) <r(s 3 ) 



+ 



s 2 2 {-8ml + 38m 2 K ) + s 2 (s 3 {19m 2 K - 4ml) 
-lUm 2 K ml - 23m 4 K + 32m^) + s 3 (16m£ - 76m 2 K )ml 
— 36m^ + \h\m\m\ + 65m%ml) a(s 2 ) 



+ (exchange Si and S2 in the second term) > , 



(D.3) 



and 



B 



(2) 

eoq-\ 



flfK 2{m\ - m\ 



{5ml -m\- 3s 3 ) , 



r (4) 



/l/x V167r(m|- - ml) 



x | _ [ S3 (8m|- - 5m2 ) + (3ml - 7m 2 K )} a(s 3 ) 

+- [2s 2 2 + s 2 (s 3 - 3(m\ + 2ml)) + (-4s 3 + 5m£ + 7m £.)] 

+i [2s? + Sl (s 3 - 3(rr& + 2ml)) + m^(-4s 3 + 5m 2 n + 7rr&] I .(D.4) 



The function a(s) appearing in all the formulas above is 



a(s) 



1 



4 m 2 



(D.5) 



In all the expressions at LO we use fx = fw = Fo- 
At NLO, we get 



\Anlo\ 2 



Bi 2) + B< 4 >) Re G 8 + (l?g> + <) G 27 + £ ^ 4) Re K t 



i=i,n 



+ 



Cg "*Re Gg + C 2 jG 2 j 



A|^ i0 | 2 = -21m G g 



^ 2) +4 4) 



G 27 ( C^? 



i=i,n 



+4?) 



(c<« + cf) 



+ fff'Re K 



-21m (e*G E ) { ( £ [Re G, (if > + Bf )] + £ ff«Re it) Cf 
I \i=8,27 i=l,ll J 

- (4 2) + 4 4) ) ( E [** <* (cf ) + cf >)] + E ^ k) ) 

\i=8,27 i=l,ll J ) 
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+ [ 2 £ ^ 4) Im £ 



i=i,n 



-2 I £ Hf } Im A'., 



. i=l,li 



j=8,27 



i=l,ll J 



i=8,27 



+ £ finite 

£=1,11 J 



Again, we disregarded the Sj dependence of the functions B^ 2n \ cf n ^ and H^ a> . The 
functions -Bg^) anc ^ ^«' 4 ^ can ^ e deduced from the results in [10] and the functions B^' 
from Subsection B.l. Finally, the functions cf^ and are discussed in Appendix E. 



(D.6) 



r(2n) 



E. Final State Interactions at NLO 

In this Appendix we provide some details of the calculation of the FSI using the optical 
theorem in the framework of CHPT. We compute the imaginary part of the amplitudes 
at 0(p e ). The calculation corresponds to the diagrams shown in Figures 1 and 2. We 






Figure 1: Relevant diagrams for the calculation of FSI for K + — » ir + ir + ir . The square vertex is 
the weak vertex and the round one is the strong vertex 




the weak vertex and the round one is the strong vertex. 

can distinguish the cases in which the weak vertex is of 0(p A ) and the strong vertex of 
order 0(p 2 ) and the inverse case in which the weak vertex is of order 0(p 2 ) and the strong 
vertex of order 0(p 4 ). In this paper we will not consider the weak vertices generated by the 
electroweak penguin. In Subsection E.l we provide some notation. In Subsections E.2 and 
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E.3 we report the calculation for the charged Kaon decays. An example of the calculation 
of the integrals that must be performed is given in Subsection E.4. Finally, in Subsection 
E.5 we give analytical results for the strong phases at NLO. 

E.l Notation 

In order to be concise we use the functions Mj for the weak amplitudes given in [10]. We 
define 

Mi{s) = J d cos 9 Mi(a{s) + b{s) cos 6)\ p4 , (E.l) 

Mf(s) = J dcos9 (a{s) + b{s) cos 9) Mi{a(s) + b(s) cos 9)\ p4 , (E.2) 

Mf s {s) = J d cos 9 (a{s) + b{s) cos 9) 2 Mi(a{s) + b{s) cos 9)\ p4 , (E.3) 
a(s) = ^(m 2 K + 3ml - s) , 



b(s) = l^j(s- 4m2) (s - 2{m\ + m£) + . ( E .4) 

The amplitudes at 0(p 4 ) for the hit — > tttt scattering in a theory with three flavors can 
be found in [59]. We decompose the amplitudes in the various cases as follows. For the 
case 7T+7T+ — > 7T+7T+ the amplitude at 0{p i ) is 

IIi = Pi(s) + P 2 (s,t) + P 2 (s,u). (E.5) 

For the case 7r°7r° — > it + it~ the amplitude at 0(p A ) is 

n 2 = P 3 (s) + P 4 {s,t) + P 4 (s,u). (E.6) 

For the case 7r + 7r~ — ► ir + ir~ the amplitude at 0(p 4 ) is 

n 3 = P 5 (s) + P 6 (s,t) + P 6 (s,u) + P 7 (s,t) - P 7 (s,u). (E.7) 

Finally the amplitude 7r°7r° — > 7r°7r° at C(p 4 ) is 

n 4 = P 8 (s) + P 8 (t)+P 8 (u). (E.8) 

The value for the various Pi can be deduced from [59]. In the following we use 

dcos# s n (c(s)(l - cos0)) m p(s,c(s)(l - cos#)) , (E.9) 

/i 
dcos# (c(s)(l - cosfl)) n p(c(s)(l + cos0),c(s)(l - cos#)) , (E.10) 



/: 



pW(a) = / dcos# (c(s)(l -cos#)) n P(c(s)(l -cos#),s), (E.ll) 



c( S ) = -±(a-4m£). (E.12) 
Another function we use in the next subsections is a(s) which was defined already in (D.5). 
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E.2 Final State Interactions for K + — ► 7r + 7r + 7r 



We first compute the contributions depicted in Figure 1 in which the weak vertex is of 
0(p 4 ) and the strong vertex of 0(p 2 ). The results for the diagrams A and B are 



Inn4 



(6,1) _ o-(s 3 ) (2ml ~ S3) 



w 



ImA 



(6,2) 
W 



32vr /2 
+3t77 2 - 2s 3 ) -M? 2 (s 3 ) 
(si -ml) 



M 10 {s 3 )\ pi + Mn(s 3 ) + Mi 2 (s 3 )(m^ 



(E.13) 



327T /2 

-M 9 s ( Sl )l + (si s 2 ) , 



M 7 { Sl )\ p4 + M 8 (ai) + M 9 ( Sl )(m^ + 3m 2 n - 2 Sl ) 



(E.14) 



respectively. For diagram C we have both S*-wave and P-wave contributions. We get for 
them 



Inn4 



(6,3) 
W,S 



aisi si 



2 M n ( Sl )| p 4 + M n (si) + Mio(si) + Mf 2 (ai) 



ImA 



(6,3) 
W,P 



64vr / 2 

-M i2 (si)(m|- + 37712 -2si) + ( Sl «->s 2 ), 
1 si(s 3 - s 2 ) 



(E.15) 



647T sf - 2(777^ + ml)si + (mj, - ml) 2 
x(M n (si) - Mio(si) + Mi 2 (si)(2si - m\ - 3ml)) + 2Mfi(«i) 
-2Mf ( Sl ) + Mf 2 ( ai )(5si - 3(7771, + 3m 2 }) + 2M«( ai ) 

+^6 2 (si) M 12 ( Sl )| p4 ] + (si ~ s 2 ) , (E.16) 

respectively. 

Secondly, we report the calculation of the case in which the strong vertex is of 0(p 4 ) 
and the weak vertex is of 0(p 2 ). With analogous notation as above, we get 



ImAl 6 ' 2 ) 



M 10 ( S3 )| p2 (Px^ + Pf ^)), 



ImA 



(6,3) 



ImA 



(6,3) 



_ cr(s 3 ) 
32vr 

= ^ (M 7 ( Sl )+M 8 ( S2 ) + M 8 ( S3 ))| p2 (P 3 ( Sl )+Pr ) ( S i)) 
+ (si ^ s 2 ) , 

= ^ (M 10 ( S3 )| p2 + M 10 ( S2 )| p2 ) (P 5 ( S i) + P 6 ( °' 0) (.i)) 

+ (si <-> s 2 ) , 
cr(si) 



32vr 



Mi (s 3 )| p2 



M 10 (s 2 ) 



p J si — Ami 



: Sl -Aml)P^ \ Sl ) 



+2P 7 (0 ' 1) ( Sl )) + ( Sl ^ S2 ). 
The final result for the Inn4( 6 ) is given by the sum 



Iim4( 6 ) 



E 

i=l,2; j=W,7r 



JmAf ' l) + 



Im . 1 



(6,3) 



(E.17) 

(E.18) 
(E.19) 

(E.20) 
(E.21) 
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The relation between this imaginary amplitude and the functions defined in Appendix 



B is 



(E.22) 



Im#= £ GiCr(si,S2,s 3 ) + Hl 6 \ Sl , S2 ,s 3 )K t . 

i=8,27 i=l,ll 

This relation is also valid for K + — > 7r°7r°7r + . 

E.3 Final State Interactions for K+ -> 7r 7r 7r+ 

The calculation is analogous to the one for K + — ► 7r + 7r + 7r~. The relevant graphs are 
depicted in Figure 2. In the case in which the weak vertex is of 0(p 4 ), we get 



ImA 



(6,1) _ 0"(S 3 ) (s 3 - 771 2 ) 



32TT /2 

+3m 2 -2s 3 ) + Mf 2 (s 3 ) 



2M n (s 3 ) + M n (s 3 ) + Mi (s 3 ) - Mi 2 (s 3 )(m^ 



ImA 



ImA 



(6,2) _ cr(s 3 ) ml 
327T /2 

-M|(s 3 ) 



(6,3) _ crQl) (2m^ - si) 



M 7 (s 3 ) + M 8 (s 3 ) + Mg(s 3 ){m 2 K + 3mJ - 2s 3 ) 



(E.23) 



(E.24) 



w,s 



2M 8 ( Sl ) + M 8 («i) + M 7 (si) - M 9 (si)(rr& 



64vr /2 

+3m 2 - 2si) +M 9 s (si)j + (ai <-► s 2 ) • (E.25) 

Also in this case diagram C generates both S'-wave and P-wave contributions. The 
P-wave contribution due to the diagram C in 2 is 

.(6,3) _ cr(si) 1 gi(g3 ~ gg) 



ImA 



[(si - (m^ + 3ml)) 



w > p 64vr f%sl-2(m 2 K + ml) Sl + {m 2 K -ml) 2 

x(M 8 (si) - M 7 (si) + M 9 (ai)(2ai - - 3m')) + 2M|( Sl ) 

-2M|( Sl ) + M|(si)(5ai - 3(m^ + 3mJ)) + 2M| s ( Sl ) 
8 

+-6 2 (si)M 9 (si)J +(si <->s 2 ). 
If the strong vertex is 0(p 4 ) and the weak vertex is order 0(p 2 ), we get 



wl 6 < 2 ) 



ImA 



<r(s 3 ) 
32vr 

32vr 

(6,3) _ g(«l) 



(M 10 (si) + M 10 ( S2 ))| p2 (P 3 (s 3 ) + Pi°' 0) (s 3 )) , 
(M 7 (s 3 ) + M 8 ( Sl ) + M 8 (s 2 ))l 2 (P 8 (s 3 ) + J P 8 (0 '° ) ( S 3)) 



; , - (2M 8 (si)+M 8 (s 2 ) + M 8 (s 3 ) + M 7 (s 2 )+M 7 (s 3 ))| p2 

x (pf°)( Sl ) + pW( ai ) + p(° 4 )( Sl )) + ( Sl «-> S2 ) , 



(E.26) 
(E.27) 
(E.28) 
(E.29) 



ImA 



(6,3) _ cr(si) 



7T,P 



64vr 



(M 7 (s 3 ) - ^s(ss) - M 7 {s 2 ) + M 8 (s 2 ))l 2 



si — Ami 

x ((si - 4m 2 )(Pr>(si) - pW (ai ) + P 2 %)) + 2P< 1 '°>(, 1 ) 
-2P 1 ( i 1 4 ) (si) + 2P 2 ( i 1 4 ) (si)) + (si ~ s 2 ) . 



(E.30) 



- 39 - 



The total contribution is given by the sum of (E.21) with the proper right-hand side terms. 



E.4 Integrals 

The integrals necessary to compute the two-bubble FSI we discussed in the previous sub- 
section can be calculated generalizing the method outlined in [60] . As an example we show 
the integration of the function 



2r„ \ „ f 2775 , rj-S A, (A - 1) 2 - tfb 2 



where Cb is a term which does not depend on t, 



C B = 2 1 - In 



r? 2 - 5 2 



and 



(E.31) 



(E.32) 



T] = mi + 7722 
S = mi — TO2 

A= Vi(t-ri 2 )(t-6 2 )}. 



(E.33) 



In the center of mass frame one can define 



Q = PK+Pn = (Va, 0,0,0) 



P „ = £ ((i - m * - ml ) , o, o, JT^ m " + gj + Kl 



m 



(E.34) 

where is the momentum of the external pion entering in the same vertex of the Kaon. 
The functions B can also be generated in the strong vertex. In this case pk is the momentum 
of an external pion. The contribution to the imaginary part of the amplitude A is 



1 f 1 
Imyl= a(s) / d cos 6 B[mi, m,2, t] 

327T J_i 



(E.35) 



with 



t = a + b cos 6 



(E.36) 



and a = a(s), b = b(s) in (E.l). In order to solve the difficult part of the integral one can 
put 



t = 



(E.37) 
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In this way 

r-l 



/- 



A (A - t) 2 - v 2 5 2 



2b 
V 2 -5 2 
2bx 



dx 



(1 — x ) lnx 



x 2 (x 2 + 1 — 2xa) 
1 — x 2 — (1 — x 2 ) In £ + ax In 2 x 



(E.38) 



-\-2x\J a 2 — 1 (in x In — 

\ I — a — 



1 — a + xV a 2 — 1 
:\/a 2 — 1 



a + Va 2 — 1 
+Li 2 (x(a + yV - 1) 



(E.39) 



where 



(V + 5 2 ) 
(r? 2 - 5 2 ) 



Xmax — 



<5 2 



if 



+ 5 2 



., - a - b + iv/(2(a + b) - ( V 2 + 5 2 )) 2 - (r/ 2 - <5 2 ) 2 } 

(E.40) 

In the case mx = rn n , a + b = and one recovers the formulas of [60]. 

E.5 Analytical Results for the Dominant FSI Phases at NLO 

The elements of the matrices defined in (6.6) have the next analytical expressions at NLO 

rLO = fRu r 12 \ 
y i?2i R22 J 

e * + ^ +) ) + e + *s )(oo+) ) * 



<5 9 NLO 



i=8,27,E 



i=l,U 



y: g, (b&+-> + B <r>) + e («. < ? ,++ -» + < 4 r +) ) * 



(E.41) 



i=8,27,E 



i=l,ll 



with 
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, NR 
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21 



xR-NR 



(/?i+/? 3 ) NR (-/3i + ^ 3 ) 
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The definitions of a±, 03, (3\ and (3s are in (6.1) and the values of their relevant combinations 



are 
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NR 
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1 
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i=8,27,£ i=l,ll 
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R-NR 



e c4r 0) + e ^? (+ - o) ^ 

=8,27,£ i=l,ll 



(E.43) 



where the functions £^0(1), Cj,0(i) an d Hi,o(i) are those obtained form the expansion in 
(B.9) of the corresponding full quantities that can be found in Appendix D. 

Disregarding the tiny CP-violating (less than 1%) and the effects of order e 2 p 2 (the 
loop contribution is less than 2%), we obtain the numbers in (6.10). 
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